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S2 Abstract 

We set up a new, flexible approach for the tensor reduction of one-loop Feynman integrals. The 5-point 
O tensor integrals up to rank R = 5 are expressed by 4-point tensor integrals of rank R—l, such that the 
^ appearance of the inverse 5-point Gram determinant is avoided. The 4-point tensor coefficients are 
represented in terms of 4-point integrals, defined in d dimensions, 4 — 2£<d<4 — 2£ + 2{R — 1), 
with higher powers of the propagators. They can be further reduced to expressions which stay free of 
the inverse 4-point Gram determinants but contain higher-dimensional 4-point integrals with only the 
first power of scalar propagators, plus 3-point tensor coefficients. A direct evaluation of the higher 
dimensional 4-point functions would avoid the appearance of inverse powers of the Gram determinants 
completely. The simplest approach, however, is to apply here dimensional recurrence relations in order 
to reduce them to the familiar 2- to 4-point functions in generic dimension d = 4 — 2£, introducing 
thereby coefficients with inverse 4-point Gram determinants up to power R for tensors of rank R. For 
small or vanishing Gram determinants - where this reduction is not applicable - we use analytic 
expansions in positive powers of the Gram determinants. Improving the convergence of the expansions 
substantially with Fade approximants we close up to the evaluation of the 4-point tensor coefficients 
for larger Gram determinants. Finally, some relations are discussed which may be useful for analytic 
simplifications of Feynman diagrams. 
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1 Introduction 



The evaluation of tensorial Feynman integrals with n external legs is an important technical ingredient 
of perturbative quantum field theoretical calculations with Feynman diagrams. They are needed in 
particular for the fast and efficient numerical evaluation of next-to-leading order contributions at high 
energy colliders. Special attention is concentrated these days on experiments performed at the LHC; 
for a snapshot on related activities see [1]. Of course, there is an unlimited variety of other reasons 
to use tensor reductions of Feynman integrals with quite diverse requirements in detail, and a unique 
all-purpose, final approach does not exist. For a recent overview see [2] and references therein. 

The first systematic approach to reduce tensor components to a basis of scalar 1 -point to 4-point 
integrals in generic dimension d = 4 — 2e for n < 4 is the Passarino-Veltman reduction [O, obtained 
by solving a system of linear equations. In fact, this is a unique basis. 

We use, with some sophistication, Davydychev's approach [4J, where n-point tensor coefficients 
are represented in terms of scalar Feynman integrals. For tensors of rank R they are defined in 
space-time dimensions up to 4 — 2e -f 2R, with an additional modification: propagators may appear 
with higher powers. These integrals are complicated objects, and an important step towards their 
evaluation is the application of dimensional recurrence relations, derived for L-loop functions in 
They have been systematically worked out for L = 1 in [6], and in a subsequent article f7] the 
evaluation of scalar integrals in d dimensions with powers 1 of the scalar propagators is advocated. 
Alternatively, the straightforward derivation of representations in the generic dimension 4 — 2£, or 
finally just in four dimensions, by means of recurrence relations introduces coefficients containing 
inverse Gram determinants, which may become small in some kinematical domains and thus raise 
numerical problems. For n <4 this problem was not very severe [8J. Serious problems arise, however, 
for n-point Feynman integrals with n>5. In that case the choice of the tensor basis is not unique 
and the freedom may be used to completely avoid the appearance of inverse Gram determinants for 
n>5||9l[IOl[ni[l2l[l3l[l4l|15l. Onthe contrary, for n < 5 one has to find explicit methods to stabilize 
the numerics for vanishing or small Gram determinants. 

In view of these facts, one may wonder if the approach of Davydychev can be worked out with an 
optimization of the handling of exceptional (small or vanishing) Gram determinants. This is what has 
been achieved in the present work. 

For non-exceptional kinematics the g^^ tensor, considering 5-point functions, is redundant and may 
be expressed in terms of 4-momenta. A nicely compact algebraic result is obtained due to this ansatz 
after applying symmetrized dimensional recurrences , but inverse Gram determinants of 5-point as well 
as of 4-point functions are introduced. In an earlier attempt [ r6l[T7l[T8l . tensor ranks until R = 3 were 
presented without inverse Gram determinants of the 5-point function using the algebra of the signed 
minors IfTQll , but a generalization to higher ranks was not evident. In the present work we first apply a 
particular recursion, which was obtained in f2\\. It reduces the tensor rank from Rto R — l and the 
further reduction can be arranged in a systematic manner without introducing inverse 5-point Gram 
determinants. The 4-point tensor coefficients, which are in fact, due to [4], higher-dimensional 4-point 
functions with higher powers of the scalar propagators, were reduced by lengthy algebraic calculations 
to higher-dimensional 4-point functions with powers 1 of the scalar propgators plus tensor coefficients 
of 3-point functions. The latter higher-dimensional 4-point functions are tensor coefficients of the 
gP-i^^i . ..g^ii-\M terms of 4-point tensors. At this stage inverse 4-point Gram determinants are avoided 
completely. The recursions of all remaining 3-point functions may be performed simply a la [6J. This 
is the way the numerics for this article was performed. Nevertheless, the same approach as above can 
also be applied to the 3-point functions, leaving only higher dimensional 3-point functions with powers 
1 of the scalar propagators, thereby avoiding inverse Gram determinants of the 3-point function. 
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Finally, one has to calculate the higher dimensional 4-point functions with powers 1 of the scalar 
propagators. This may be done by direct evaluation, see e.g. [|22|. or by further reduction to simpler 
integrals, which in general, however, introduces inverse powers of 4-point Gram determinants. For 
small Gram determinants we therefore derive a relatively simple analytic expansion in positive powers 
of the Gram determinant. Such an infinite series was, to our knowledge, first proposed in equation (36) 
of [Ell, but was not numerically applied so far. In fact, this expansion applies for higher-dimensional 
integrals with powers 1 of the scalar propagators only and would not be appropriate for a representation 
of tensor coefficients in their original form. 

Another approach to the problem of small Gram determinants was chosen in sect. 5.4 of ['231, 
where relations between different 3- and 4-point tensor coefficients are exploited and the full set is 
calculated with increasing iterations. The n''^ iteration requires all 3 -point coefficients of rank n. In 
contrary to that our expansion is concerned only with the subset of the g^^^^ ■ .. gl^ii-\^n coefficients of 
the 4-point tensors, which are approximated only by the corresponding subset of the 3-point tensor 
coefficients. Indeed, our series expansion of the higher dimensional 4-point functions is useful only 
since these are embedded in expressions which are already free of inverse Gram determinants . 

In [[7]| it was shown that the integrals under consideration can be expressed in terms of multiple 
hypergeometric functions. One can apply their series expansion or, alternatively, one could use their 
representations in terms of 1 -dimensional integrals also given in [[Tj. These, in general, present the 
integrals in rather different domains of phase space, including the case of small Gram determinants of 
the 4-point function. Thus our approach offers a variety of options to adjust to the given kinematical 
situation. For the time being we only use our series expansion applying Pade approximants. This 
turns out to be very efficient and allows to obtain high precision for the numerical values of the 
tensor coefficients of the 4-point function. In an example, we demonstrate that the combination of 
representations for non-exceptional kinematics with this expansion covers the complete phase-space 
from medium to vanishing Gram determinants. 

In view of the importance of stable numerics for tensor reductions, it would be welcome to have 
one or more complete opensource programs for this task, including the treatment of small Gram 
determinants. To our knowledge, none is presently available. Following the approach of this article, a 
C-i-i- program is under development to close this gap ['241. 

The article is organized as follows. In sect. [2} some definitions and basic formulae are recalled. 
Sect. [3] describes a compact analytical tensor reduction of 5-point functions with non-exceptional 
kinematics. In sect.|4]the 5-point functions up to ranki? = 5 are reduced to 4-point tensor coefficients 
in terms of 4-point integrals in higher dimensions and with higher powers of the scalar propagators. 
In sect. |5| we reduce these 4-point integrals in several steps. In sect. 5.2 the integrals are reduced to 
4-point integrals in higher dimensions with powers 1 of the scalar propagators plus 3-point tensor 
coefficients. The results are given in eqns. ( [5.10[ ), ( [5.15[ ), ( [5.19[ ) and ( |5.21[ ). Indeed, these eqns. are the 
central point of our approach since they allow to proceed further in different directions. One might e.g. 
apply the general method of calculating higher dimensional integrals of [7 [. With app.|Aj alternatively 
a reduction to the Passarino-Veltman basis is straightforward. In subsect. 5.3 we recall how to expand 
the higher-dimensional integrals with powers 1 of the scalar propagators for the case of vanishing 
or small Gram determinants, the result being given in (5.46). The symmetrized recursion relations, 
useful for the 5-point functions with non-exceptional Gram determinants are presented in sect. |6} 
Additionally, in sect. |7] we give some relations which may be useful for an analytic simplification of 
Feynman diagrams. We end with conclusions in sect. [8] Appendix |A| contains a list of dimensional 
recurrences and app.|B] collects divergent parts of higher-dimensional integrals. A numerical example 
is discussed in app.|C] Appendix |D] contains notations and some relevant algebraic relations. 
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Figure 2. 1 : Momenta flow of the «-point function. 

2 Tensor integrals in terms of integrals in shifted dimensions 



A tensorial Feynman integral with n external legs is shown in fig. 2.1 and is defined as 



n,{Vj} ij^2 

with denominators Cj, having indices Vj and chords qj 



J = {k-q.f-m] + ie. (2.2) 

Here, we use the generic dimension d = A — 2e and /i = 1. Reducing the tensors to 1- to 4-point scalar 
functions in general their expansions in terms of e is needed. The first expansion terms can be 
expressed in terms of Euler dilogarithmic (or simpler) functions [|26ll27ll28ll25ll29l . 

The six-point tensor integrals may be expressed in terms of five-point tensor functions [TTl |30l 
21: 
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where the auxiliary vectors Qs are 
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i=i We 



A similar formula exists also for five-point tensor integrals [|2T]| : 



The auxiliary vectors here are: 



Qs = i^flr' ^ = 0--'5. (2.6) 



For later use, we introduce also 



/.A 

uVah 



(2.7) 



a 



tu,il 



5 

EJi\itJ_5 
Hi nu\ 
1=1 \tu)5 



(su\ 
\tu) 5 



\tu)' 



and 



5 ('■) 



:i ()5 



(2.8) 



(2.9) 



In fact, ( |2.3[ ) is essentially the same formula as ( |2.5| ), except that ()g is replaced by (q)^ etc. and 
(00) = 0. With the definition 



m, 



(2.10) 



the modified Cayley determinant of a topology with internal lines 1 • • • n becomes 



0. = 



1 1 ... 1 

1 Yi2 ■■■ Yin 
1 Yn Y22 ... Y2n 



(2.11) 



1 Yin Y2n ... in- 
One chord may be chosen arbitrarily to vanish, qn = 0, and then this object is the Gram determinant: 



On\q„=o = — det G„_l, 
Gn-i,ik = '^qiqk, i,k=\,...n-l. 



(2.12) 
(2.13) 



The Gram determinant is independent of the internal masses. The signed minors [fT9ll are denoted as 
follows: 



Jl ]2 ■■■Jm 
ki ^2 ■ ■ ■ km 



(2.14) 



They are determinants, labeled by those rows ji , j2, ■ " ■ jm and columns ki,k2, - ■ - km which have been 
excluded from the definition of the Gram determinant ()„, with sign 



Jl J2 ■•■ Jm 

m 



sign(- - -J =(-l)^'+^'^+---+^"'+''+'^+'''+'"-Sai,72---jg-5(A:i,A:2r (2.15) 



'Usually we will use indices 5,f , • • • = 1, • • • ,« for labelling internal lines, and indices i,],-- - — I.-- - ,n—lfoi labelling 
the (non- vanishing) chords. 
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Here S{ji , j2 • • ■ jm) gives the sign of permutations needed to place the indices in increasing order. |^ 
We have e.g. 



A„ = 



Yn 


Yn . 


■ Yin 


Yu 


Y22 ■ 


■ Yin 


Ym 


Yin . 


Ynn 



(2.16) 



Applying Davydychev's method [HI, one expresses the tensor integrals ijf^"'^" by scalar Feynman 

integrals ijf^,,, in higher dimensions d and with higher indices v,. We reproduce here integrals with 
ranki? < 5: 



vA 



^.^^ n ^7 



-1 



7=1 
nA ' 



(2.17) 



n 1 



(2.18) 



L ^j^k'^ijk^n,ijk +jLS ^n,i 



d'^k 



(2.19) 



2^ 'ii 'ij^k'il ^ijkl ^ndjkl 2^ 1i1j^iJ^n,ij S 



ij,k,l=l 

d'^k 



(2.20) 



E" M V A P (7 Ad+f s ^ \^ [ixv A P C7] Ad+f 

Qj Gk nijklm In,,jklm +3 ^ ^ ^' "'Mvji 

',j,k,l,m=l ij,k=l 

-^^;g[^v^^^^^^^ (2.21) 



^The definitions are related to similar ones used in the literature, see app. |d] 
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The following symmetrized tensors are used: 

= g^'qf+g^Ul+g'^ql (2.22) 

g^^'qUf = g^'qU'j+g''^qU'j+g'^qU' + s^'<iU]+g''qU]+s^'^U]^ 

(2.23) 

^[Mv^Ap] ^ (2.24) 

g^'^'qWjqf = g''''qfq'jq!+g'''qJq'jq'k+g''qrqjqk+g'''qrqjq1+s'"q^q'jq^ 

+ g''qU]4+g'''qU]€+s^'<iU]4+g^''qUK+s'''^UlAA2.25) 



The scalar integrals are: 



[d+tstu- ^ / ■ - T-f t (2 27) 

T — 1 Cf 

where [<i+]' = 4 — 2e + 2/. The index p is the number of propagators of the p-point function. Note 
that equal lower and upper indices cancel. The coefficients nij.riijk and riijti etc. in (2.18) to (2.21 ) 
were introduced in [TTtII . They stand for the product of factorials of the number of equal indices: e.g. 
nun = ^^-1^1 jii = 3\,niijj — 2\2\, niji^i^ = 2\,niji^i = 1 !; the indices i,j,k,l are assumed here to be different 
from each other. The following relations are of particular relevance for the successive application of 
recurrence relations to reduce higher-dimensional integrals: 

nijk = ^ij'^ijki 

nijki = "^ij^ijk^ijkh 

nijklm = '^ij'^ijk^ijkl^ijklmi (2.28) 

and: 

^ijk = 1 + Sik + djk, 

yijki = '^ + Sii + dji + dki. 

^ijklm = ^ + Sim + Sjm + Skm + Sim- (2.29) 

In a second step, one may choose to express the higher-dimensional scalar integrals in terms of 
the generic scalar integrals The algorithm is based on recurrence relations with shifts of dimension 
J > 4 — 2£ and indices v.^ > 1 , 

(),v,(s+/,r^') = +E(;)jt-/<'"), (2.30) 

or with a shift of dimension d: 

()„(d-tv,+ l)ri«> = Qf-t(°)^^{t- !'/>)■ (2.31) 



These relations hold for arbitrary index sets {v,}. The integrals s+/„ and t /„ are obtained from 
^^by: 
app.g- 



ijf' by replacing v., — )■ (v^ + 1) and — t- (v^ — 1), respectively. For more explicit expressions see 



3 An efficient reduction of 5 -point tensor integrals 



The reduction of 5-point tensor integrals to 4-point tensor integrals at non-exceptional momenta may 



be performed by iterative application of ( |2.5[ ). This was exemplified in pT|, and an opensource Fortran 
code olotic ||3T1l is available. In this sect., we derive a very compact, explicit representation of the 
tensor coefficients for 5-point functions in a minimal basis, chosen to be free of the metric tensor. This 
will rely on an exploitation of ( |2.9[ ) and ( |3.9[ ), a specifically useful relation of the algebra of the signed 
minors, and applying the symmetrized dimensional recurrences of sect. |6] 

We investigate the 5-point tensor integrals step by step. For the tensor of rank i? = 1 we get 
from (1231) 



5=1 



(3.1) 



and /5 may be taken from ( |4.2| ). 

Similarly for the tensor of rank 2, 



(3.2) 



with four-point integrals from (2.17) and (A. 6), 



^4 



such that we can write the tensor of rank R 



\is/5 ts I V \is/5 jst 

= 2 as 



,t-, 



3 : 



s)5 



(3.3) 
(3.4) 



^5 



t=l 



(3.5) 



Compared to ( |2.18[ ), this representation and the following ones are free of the metric tensor. Further, 
the compactness relies on the use of the auxiliary vectors Qs ,Qt '^ instead of the chords . 
The tensor of rank R = 3 deserves a bit more effort. 



(3.6) 



The corresponding 4-point function reads now due to ( |2.19[ ) and with ( |A.5[ ) 



^4 



W4,ij 



1 



_ IXVT[d+],s 
2^ 4 ' 

_\jsl5_j[d+],s \jsJj_ [d+]j ^ \jsJ_5j[d+],st 



(3.7) 



(3.8) 



Observe that for j = s the integrals I^^ and -j ' vanish (due to vanishing signed minors) such 



that indeed a formal summation over all five values of ?, j is possible. 
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Now we use identity (|D.3|), 



Os O5 ()5(;:)5 

()5 Os ' 



(3.9) 



where the Ql are the vector components of Qs, see ( |2.6| ). 

Performing summation over ij in ( |3.7[ ), the first term on the right hand side of p.9[ ) yields 

^g'^^/j^^''^ (see ( |2.9[ )) and thus cancels against the last term in (3.7). Thus we can write the 4-point 
tensor of rank i? = 2 as 



T^^v,s 

^4 



■ld+],s 





ijf}_5j[d+],s _ Osij^s j[d+ls ^ ijs)s j[d+],st 
■■ M + ^s^ ^3,i ' 



(3.10) 
(3.11) 



'5 O5O5 t=\ \s)5 

where the metric tensor has again disappeared compared to ( 2.18[ ) and instead q^q^j contribute for 
iJ = s. A compact notation can now be used with ( |2.7[ ) and ( |2.8| ): 

5 



4 -Qo ^0 k - J^QsQsU +hli'ij^3,K 



(3.12) 



7=1 



withi?!; ' the scratched version of (6.2). Inserting (3.12) in (3.6) yields the compact expression for 



,[d+f,s 
'3,0- 

I^^^, free of the metric tensor. 



Next, for the tensor of rank 7? = 4 of the 5-point function, 

5 



^5 ~h '^0 Z^h 



P 



(3.13) 



.5=1 



we need the 4-point function of rank = 3, according to ( |2.19| ): 

5 1 5 



(3.14) 



i,j,k=\ 

lQls\ lOs\ tQs\ 
\is) 5\is' 'i.\ks) 



•'+u"*)+('"*) (3-15) 



0; 







sJ5 



Here we have used for 7^ -^^ '' , instead of ( |A.4 ), the symmetrized form ( |6.8 ) with i?^ -^^ ' being the 
scratched version of \6.9) . Applying again p.9[ ) we obtain the analogue of ( |3.12[ ), 



•'A,ijk 



I ^UhkJl,ijk. 

{is)5ijs) ^ik^5 1 



(3.16) 



0; 



0^ 
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' + {j^k) + {i ^ k) } -<+i '\ (3.17) 



3,ijk 



For the following it also pays to introduce 



t=i 



so that finally 



5 



(3.18) 



(3.19) 



with ^1^,^], the scratched version of (|6.9|). This finishes the determination of I^'^^^ . 



Finally, for the tensor of rank i? = 5 of the 5-point function, 

5 

E 

^=1 



(3.20) 



we need the scratched tensor of rank i? = 4 of the 4-point function. The corresponding symmetrized 

tensor coefficients are taken from ( |6.11| ) and ( |6.1 1[ ) by scratching. We begin with the term I^^^ 
from ( |6.11| ). Using again p.9[ ) we have 



{ls)5 (Ds I (is) 5 (fa) 5 I (/^)5 O'^-^S ^ (/)5 (^5 , (Ds (^5 , (/)5 t 



Os (::)5 Os Os Os Os Os Os Os Os Os Ot 

Os ) ^illnJnk^^n'n' 



Os I Os 

Os 



\>s y>s Us Us Us 



(3.21) 



The first term on the right hand side of (3.21 ) yields after summation over /, j, k, I 



Lfi V X P 



(;)sa)s , (OsQs , (t)sC)sl_i 



+ 



+ 



ijkl=\ 



Os Os Os Os Os 



(3.22) 



The same contribution comes directly from ( |2.20| ). Finally, the second term of ( |2.20| ) contributes with 

j[d+]\s 
^4 



the term ~ /J, of (6.11 1 and again the first term on the right hand side of ( 3.9 1 



2 £^ o Hi Hj Q ^4 ^6 6 ^4 ; 



E 



(3.23) 



i.e. the terms g^^^ g 



[pVpAp] . ^^jj^gi j]^g second term on the right hand side of p.21| ) yields after 
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summation over /, j 



V5/5 ij.k,l=l 



(3.24) 



term^/J^+J'''of([6ll]), 



A contribution of this type also comes from the second term of p.20[ ) with the 
but now the second part on the right hand side of p.9| ) is 

which means that also these terms cancel. Finally there remains the last term in ( |3.21[ ), which does not 
cancel but contains no g^^ 



(3.25) 



Us 



(3.26) 



The contributions from ( |6.1 1| ) of the type jff^ - after cancelling the first term on the right hand side 
of ( |6.11| ) - can also be written in a compact manner, 



IJ5 



(3.27) 



where the symmetrization in the tensor indices is understood. Introducing like in p.l8[ ) 



,stu 



(3.28) 



we can finally write 



Os Us 



t=l { [stJs 



i,j,k=l t=l ij,k=l 



(3.29) 
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They;-' defined in (3.18) occurs again, the Rl I 



is given in and Rj'tp" in ( [6l3] ). With the 



expression for free of the metric tensor, we complete the rewriting of i^^^P^ with p.20[). 



It is remarkable that all the coefficients of the g^^ terms of the four-point functions in ( |2.18[ )-( |2!20] ) 
completely cancel in a way that the remaining tensor coefficients are much simpler than the original 
ones. This is achieved due to the symmetrization of the recurrence relations given in sect. [6] and would 
have been seen less easily with the "standard" recursions of app. |Aj The new representations for the 
tensors may be useful in several respects. First of all we have here an extremely compact notation, due 
to the use of auxiliary vectors Qs, which is not evident at the outset. Further, the representations may 
be used for a completely independent programming and thus for stringent numerical cross checks. The 
latter one is an important aspect because there are not too many opportunities for that in case of the 
5-point and 6-point functions. Finally, the auxiliary vectors have some specific properties so that 
they may be used for simplifying manipulations with physical amplitudes, see sect. |7J 



4 Reduction of 5-point tensor coefficients 

The purpose of this sect, is to express the 5-point tensor coefficients in terms of 4-point tensor 
coefficients, which will be evaluated in Sec. |5]in such a way that also the case of inverse sub-Gram 
determinants can be dealt with in an elegant manner. The difference to the former sect, is the fact 
that we avoid all inverse Gram determinants, 1 / ()g as well as 1/ ()4. In this case we have to keep the 
^'^^ -terms. 

4.1 Scalar and vector integrals 



For the scalar 5-point function I^, we use the recurrence relation ( |2.31[ ): 



c-myr-Qj^-LQji (4.1) 



The integral I^^'^^ is finite for d = 4, and we get in this limit: 

1 ^ /O 



(0)5 -^=1 ^^^5 

i.e. the scalar 5-point function is expressed in the limit J — )• 4 in terms of scalar 4-point functions, 
which are obtained by scratching in the five terms of the sum the s''^ scalar propagator, respectively. 
This was already derived in [fT9l , see Eq. (6.1) the re. Se e also [|20ll . 



obtain quite similarly 



The tensor w-point integral of rank i? = 1 in (12. 171) can be expressed by integrals and we 



= (4.3) 

= (d+i-„)|./K+i ' (4.4) 
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where again for n = 5 in the limit J — )■ 4 the scalar integral I^^^^ disappears: 

Ei = t^l 

/OA 

pS \0sJ5 TS 

4.2 The recursion formulae 



(4.5) 
(4.6) 



For the general case, we use as a starting point ( |2.5[ ). In order to solve this eqn. recursively, we multiply 
it with (q) g l^and use the identity 



()5 + 



(4.7) 



The first term on the right-hand-side can cancel already a Gram determinant and the second one 
transforms a vector Qs into a vector Qq. Asa result, we get from ( |2.5[ ) the general form 



J 



(4.8) 



with: 



j'Mi ■••Mfi-i 



jfll...^R-US 



(4.9) 



and 



5 

;=1 



,5. 



(4.10) 



The barred vectors are free of the inverse Gram determinant Q^. Evidently, the reduction of j^^'"^'^-'^'^ 
is also free of ()5, and we have to care only about the product rMr- Ms-ig^ 

The following observation will prove to be useful: r'^i •■^^ i contains general tensor structures as 
given in ( |2.17| )- p.21[ ) with chords qt and the metric tensor. In fact, when calculating the 5 -point tensor 
recursively, we keep at this stage the 4-point tensor as given there. With the 5-point tensor of rank 
R= I, given in ( |4.4[ ) above, the recursion is started. 

In order to cancel l/()5, in each recursive step a term (j)^ will be generated and summed over with 
the corresponding chord qt. We will apply to such terms the identity 



()5 



+ 



()5 



(4.11) 



The ratio (y)g/()5 comes from Qq, see ( |2.6[ ) . In the first term of the right-hand-side of (4.1 1 1 the Gram 



determinant has cancelled and the second term yields a contribution according to (|2.9[). The 



metric tensors in the original r^^i remain unchanged. From the following examples the scheme 
will become more evident. 



of 161). 



^Throughout the present work we assume (JJ) ^ ^ 0. In case of vanishing and/or small (q)^ see the discussion in Sec. 2.2. 
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4.3 The tensor integral of rank R = 2 



Equation (4.8 1 reads for the tensor of rank R = 2: 



^5 



.?=1 



(4.12) 



.5=1 



The square bracket, a special case of ( |4.9[ ) for R = 2, will be rewritten now. We use ( |2.17[ ) and ( |4.6[ ) for 
isert 



d — 4 and insert the reduction ( |A.6| ) with I — I: 



(4.13) 



5 

i=i 

5 
i=l 



Using further 



5 



j[d+ls 
4,i 



^4 + 



Z5 



we see the cancellation of E^. Additionally, it is 



5 f=l,f/i' 
05 



jst 
^3 



(4.14) 



(4.15) 



(4.16) 



Here the (^) ^ term cancels and the remaining factor (q^) ^ is antisymmetric in 5, t, yielding a vanishing 



contribution after summation over s, t. With ( A.7 1, reintroducing , we obtain 



!=1 



i •> 



(4.17) 



(4.18) 



Here we observe the first occurrence of a term Q^, as mentioned in sect. 4.2 Using (4.11 ) and the 
notation 

5 

(4.19) 



we finally get, taking into account ( |2.17[ ) for n = 4, the expressions for the tensor coefficients: 

Eqo 



s 

00 



^=1 



3 1 1 

^ 27o^~ 



■\d+],s 



(4.20) 



s=l 
5 

E 



s 

'j 

s=l 
5 1 



^=1 (0)5 



Os 
Oj 



j[d+],S 



(4.21) 
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The functions /|'^+]''' /|'^+]''' ^[\\ ]jq further treated in sect. 



A comparison shows that the tensor 



coefficients Eij given in Eqs. (3.10)-(3.12) and (A. 22) of fut are much more involved 



4.4 Reduction of integrals with rank R = 3 



For the tensor integral of rank R — 3, Eq. (4.8 1 reads 



0\ ^vX 
0)P 



(4.22) 



Investigating the square bracket, i.e. the tensor (4.9 1 for i? = 3, we see that the corresponding g^^ term 



vanishes. Indeed, from (|4.20|) and (|2.18|) we have: 

1 IS 



^00 + 



0. 



(4.23) 



This is interesting in view of our general scheme, which was described in Sec. 4.2: Since there is no 



vector qj in this contribution, no (^)g is produced, and if we assume that no inverse Gram ()g should 
occur in this case, the contribution must vanish. 



Further, from (4.21 1, (2.18) and (A.5) we obtain 



i ^^^^^^ 



(4.24) 
(4.25) 



and 



Oi 



V5 iVJJs 
s 




0^ 
0/ 



Os 



1,.^ / IS 



MA 



With (4.15) and 



S\ I IS 

W 5 



J'/ 5 



Os 



ts 



5 \py 5 



+ 



t=l,t^s,i \J^/ 5 



s\ fOi 
W/s 



[d+],st 







(4.26) 



sJ5 



(4.27) 



we see that the complete term E-j cancels. As above we use again (|4.16|) and with the same arguments 



as before we see that only (j)g-type terms remain such that (4.11 ) can be used again to cancel the 
Gram determinant. 

Before collecting all contributions, we would like to point out that, after the above manipulations, 
the expressions are in general not explicitly symmetric in their indices, although the original integral is 
symmetric in /i, v, A. Consequently, our result must also be symmetric in the indices /, however. 



after summation over s and t. For an explicit example see also the discussion after (6.1 ). If there is no 
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explicit symmetry before summation over s and t it may be useful to symmetrize the result. With this 
in mind, collecting all contributions, we have 



'J 



+ 



05 

5 

E 



+ 



ts 



Os 

s 



■\d+],s 



5 ViV 5 



+ 



ts 

J J 5 V'V 5. 



-2 







(4.28) 



To obtain this result, the vector integral 7^, represented by tensor coefficients Z^",^'"', and the vector 



sJ5 



[d+U 



integral I^, represented by tensor coefficients I^P'^' (|4.26|), have been reduced to scalar 2-,3-, and 
4-point integrals in generic dimension d by means of (A3) and (A. 10). The 2-point functions cancel 
here. Further we need the identity 

2 

(4.29) 



Jd+],st . 



Ost 
Ost 



and in order to get rid of the vector indices in the 2-point functions, we need the relation 



ust 
jst 



ts 



ust 
Ost 



us' 

5 VJV 5 



ts 



(4.30) 



which shows that after cancellation of {f^)^, Eq. (4.30) is antisymmetric in t and u such that it can be 

effectively considered to vanish after summation over t and u. This allows finally to introduce l^'^^'^' 
according to (|A.11[) into (|4.28|). 



There is a further subtlety concerning (4.28). The ultraviolet (UV) divergency 

1 

'2£' 



jld+],St 

'3,UV 



(4.31) 



when combined with = 1 — e, yields a constant finite contribution ir Since, however. 



5 



1=1 



ts 



0, 



this term does not contribute and we can put J = 4. In that case (4.28 ) reads 



I J 



+ 



J/ 5 



r[d+r,s 



(4.32) 



(4.33) 



to be compared with (4.18). According to our general scheme, each ^, generates a factor (^)g, the 
further factor being a higher-dimensional integral with index (indices) being the same as in the 
remaining chords. In fact, (4.18 ) is a vector coefficient so that no additional index is available and thus 
the higher-dimensional integral c annot carry an in dex. 

We just mention that, due to (4.32) and (A. 6), also the integral /F^^ is UV and infrared (IR-) 

finite. Applying (4.11) in {4.22), we obtain products (0)5/4, , for which we can write, using 
( |4.15| ), ( |4!l6| ), ( |4.27J and (|A.7fand setting d^d + 2: 

'Os^ 



j[d+]\s 



0/ 



r[d+U 



{d-l)ll 



(4.34) 



^See the discussion after ( A. 1 1 
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Collecting all the contributions, our final result for the tensor of rank R = 3 can be written as follows: 



5 

I 



5 

E 

k=l 



00k, 



with 



^00; = 



E, 



ijk 



^=1 (0)5 
5 

5 1 
" E 70^~ 



1/0^ 
2V0j 



d—\fs 



i/5 



O7 



+ 



Ok 



In ( |4.36| ), we can put d = 4 because, similar to the discussion by means of (4.32 1, it is 

= 0. 



5 

E 



(4.35) 



(4.36) 



(4.37) 



(4.38) 



Another possibility to argue uses that, due to (A. 7 ),l — 2, the d — \ cancels. 

The rank R = 3 tensors were also treated in [17] . We proved there successfully the cancellation of 
1 / ()5, although the corresponding formulae were quite a bit longer than here: see Eqs. (3.41)-(3.42), 
(3.30)-(3.33), (3.40) in [17J. For the ranki? > 3, however, the tensor reduction would become really 
awkward with the older approach. 



4.5 Reduction of integrals with rank R = 4 



For the tensor integral of rank R = 4, Eq. (4.8 ) reads: 



(4.39) 



Here I^^^ is given in (4.35) to (4.37), 7^^'^''' in (2.19), taken at n = 4. In a similar manner we 
decompose the square bracket in (|4.39): 



according to which: 



^5 



^OOi 



s=l 

^j^k^ ijk + l^S Qi OOP 



5 

E 

1=1 



^00/ 



1 fs 



2 \0 



2\i 



(4.40) 
(4.41) 



(4.42) 
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Obviously, the tensor coefficient Eqqj^ has been completely eliminated - as observed before in ( |4.14| ). 
As in (|4.35|), we now write 



I q'^q';qiq?E,jkl+ I, 8^^'qhfEoo.j+g^^'g''%000. (4.43) 



i,j,kJ=l 



Proceeding as before, using ( |4.11[ ), from ( |4.42p we obtain: 

5 

£'0000 = 

1(0) 



0000' 



(4.44) 



^0000 



4 /0\ M 



5 r[d+f,s 



(4.45) 



We remark that since (4.42 1 is summed over s, all (constant UV-) divergent contributions from 
can be dropped; see a 



SO the discussion at the end of sect. 4.4 
In the next step we calculate T^j^: 



+ 



Q 1 ^ij^ijkl^Sjk 



sk , 



Os 
Ok 



IS 



ksj^^^j 



ks 



+1 



ts 



Os 

ksJ^^'J 



+ 



5 

E 

5?=1 



'd+]^S 



ks 



j[d+]\st 



□ 5' 



(4.46) 



■\d+]\s 



where again (4. 15 ), (4. 16 1 and d4.27 ) have been applied and Vjjkl^ ^jl '' has been replaced by means 
of (A. 4 1. Again we observe that the complete tensor of lower rank (here E-jj^) cancels. After further 
lengthy manipulations and subsequent symmetrization, the following analogue of ( |4.18p and ( |4.33p 
can be verified: 



Ujk 



^jk\jk + 



i/5 



^ikL 



[d^ 



A.ik 



+ 



(4.47) 



Using again ( 4.1 1[ ), we can immediately write down the pure spatial components: 

Eijki 



5 

5 1 

E 



.5=1 



Ok 
si 



'd+Y,s 



+ {i ^k) + (j ^ k) 



+ 



_^ijk^A,ijk 



(4.48) 



For the mixed terms Eqq - ■, i.e. those containing the metric tensor, we have contributions from different 
origins. From TL, (see (|4.42l)) we get 



1 



2 E S'^^'^k^i'i V ,/ 



Ok 



(4.49) 
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From (|4.47|), we get 



[d+]\s 
4Jj 



1 fs 



2 VO 



(4.50) 



Finally, there is a contribution from the second term of ( 2. 19 ): 

5 



Collecting these contributions without symmetrization we have: 



(4.51) 



^000' 



EE. 

s=l 



000' 



OA [d+]\s (Qs\ \d+]\s (S\ [d+]\s 



(4.52) 
(4.53) 



A general comment is in order at this place: The only UV divergent term in (4.53 1 is ^ ' , which 



■\d+]\s 



comes from ( |4.42| ). We see, however, that due to ( |4.38[ ) this (constant) term does not contribute when 



Tqqj^ is summed over s. Thus, the UV divergent part can be dropped in if^^ and as a consequence it 



does also not appear in (4.53 1. This, after all, is only an expression of the fact that the original tensor 



integral under consideration is finite. 



4.6 Reduction of integrals with rank R = 5 



For the tensor integral of rank R = 5, Eq. ( |4.8[ ) reads: 

5 



0\ jUvXpo 



0\ jlLvXp 



g„ _ ^ ,,vAp,.gO.<, (4 54) 

S=l 



Writing the square bracket in ( |4.54[ ) as 

5 

E 



rpllVXp _ V-l rplJ.VXp,S 



5 5 

^5 - 1^ 1i Ijlk^l ^ijkl+ 2-1 S ^iGj^OOij+S 8 -'0000' 

i,j,k,l=\ i,j=\ 



(4.55) 
(4.56) 



we, first of all, observe that the gl'^^g'^P] term vanishes. Indeed, from (4.45 1 and (2.21 ) we have: 

'0 



7i)W-(oj^^oooo-4(oy^M 



0. 



(4.57) 
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Again we have a situation like in (4.23): There is no vector and no ('1)^ is produced. Thus, no 
inverse Gram determinant appears smce this term vanishes. The next term, Tqq-j, is calculated similarly 
as is scetched in sect. |4.4| with a result generalizing (|4.42[): 







1 fs 



8 V/ 



JJs 



(4.58) 



In contrary to the discussion at the end of sect. 4.5 , summing ( 4.58[ ) over s, the UV divergence of the 



integrals I^^^^ does not drop out since I^p = for 5 = i. The corresponding divergence cancels in 
this case against a divergence coming from the last term of ( 2.20| ). 



In analogy to (4.18 1, (4.33 1 and (4.47 1 we also have 



_^jkiU,jki + 



■\d+]\s 



i/5 



j[d+\\s 
I j ^i0A,ijk ■ 



(4.59) 



It is interesting to note that a second chain of tensor coefficients has developed for the square bracket 
tensor r'^i ••^^ i (4.9) which follows the same rule when proceeding to higher ranks, namely (4.42) 
and ( |4.58| ) to be compared with the chain ( |4.18| ), ( |4.33[ ), ( |4.47[ ) and ( |4.59[ ). 
The complete tensor of rank R = 5 ( |2.21[ ) now reads 



(4.60) 



5=1 



5 5,5, 
i,j,k,l,m=l ij,k=l i=l 



(4.61) 



Using (4.10) and (4.1 1 ), we obtain for the pure spatial part 

1 rr/oz 



^ijklm 



©5 



+ 



Next we consider again the mixed terms and begin with ^ooiji- From (4.58 ) we have 

5 



From (4.59 ) we get: 

5 



(li (Ij (Ik (ll 1m \f.j \ -7^^jkl\jkl + -T^^iklU.l 
i,j,k,l.m=\ V^/S 1^ Us Us 



OA , /OA uM\s 



sk 



(4.62) 



(4.63) 



{m)t 



im)'. 



j[d+]\s 
■Jkl 



\m)t 



■ \d+]\s 



iin)' 



[d+]\s 



oM^ 1- "ijikii^jkAJi +g L qUkiinikih 



1 ^5 

2 



AM 



I A(7 V P , oa \^ U V X T\d+ 



i,j,k=\ 



[d+\\s 
ijk 



(4.64) 
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There is a 4-point contribution from the second term of (|2.20[): 



i,j,k=l 



Ok 



(4.65) 



In (4.63 ) and (4.65 1 we have the tensor structure g^^^q^qjq^, and in (4.61 1 the structure g^^'*' qf q'^j q^ • 
In order to identify them, we can make use of the fact that the tensor /^^ ''^ is symmetric in all indices 



due to which we have only to count the number of terms in the structures to be compared: e.g. (2.25 



contains ten terms, while (2.23 1 contains six terms. Thus, replacing (2.23 1, multiplied by q^, by (2.25 



we have to introduce a factor taking care of the ratio of the numbers of terms in each of them. Similarly 



this applies for (4.64 1. In this way we obtain 



^OOijk 



1 



20. 



d /Oi 



[d+]\s dfOj 



4A 



Os 
Ok 



+ 



(4.66) 



This concludes the ^oo/;! and finally we have to collect the contributions to ^oooo;- They come from 
the last term of ( |2.20| ) and ( |4.58| ) with the result 



^0000/ 



+ d 



r[d+]\s 

M.I- 



(4.67) 



We will not proceed here further, but by now it may be evident to the reader how to treat tensors of 
higher rank. One remark, however, is in order: identifying the if^} in ( |2.17[ ) - ( |2.21[ ) as 4-point tensor 



coefficients, the above tensor coefficients of the 5-point functions should be equivalent to (6.17) -(6.21) 
of [23J. Nevertheless, working with higher dimensional 4-point functions and in particular using the 
algebra of the signed minors appears advantageous to us. 



5 Calculation of higher-dimensional 4-point functions 

In the foregoing sects, the 5-point tensor coefficients . ' have been rewritten in terms of 4-point 
tensor coefficients. The factor 1 / ()5 has been completely avoided. In detail we have: 

• Tensors with R = 2: 

The tensor coefficients Eoo,Eij axe expressed by 

• Tensors with i? = 3: 222 

The tensor coefficients £00^5 -E/jyt are expressed by /j^^^''^/]^^' '^^^ifj^^ '^14'/^^ 

• Tensors with R = 4: 2233 
The tensor coefficients £0000, £'oo/y,£,7« are expressed by /j'^^' '^4^7^ ''''4^ 

• Tensors with 7? = 5: 2334 
The tensor coefficients £oooOn£'ooi;yt,£'o-«m are expressed by I^^^^ ''\ /j^^^ I^f^^l 

r[d+]\s 
\ijkl ■ 
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It is our goal to find a representation of these integrals which is suited for the most problematic cases 
occurring in practical calculations, namely for vanishing sub-Gram determinants ()4 = In the 
numerics we will make use of opensource programs for the calculation of few master integrals, chosen 
here to be the scalar 1 -point to 4-point functions in generic dimension d = 4 — £,'m standard notation 
the integrals Aq, Bq, Cq and Dq. They are available from e.g. the LoopTools/FF package [I25ll28l or 
from the QCDloop/FF package [|29ll28l . For this purpose, we have to reduce dimension and indices 



of the above integrals. This may be done by recurrence relations ( |2.30[ ) and ( |2.31[ ), given in detail in 
app.|A]. In each recursion step an inverse power of ()4 is generated, which causes numerical problems 

for small ()4 although the original integrals I^^^} are finite and well-behaved there. 

We proceed in two steps. In subsect.|5T[ an intermediate step, we manage to write the integrals in 
the form: 



[d+]',s 



\ks)5 



[d- 



,s _ [d+ 



(5.1) 



,[d 



,4 is constructed such that in the limit ()4 — )• it has the same value as 



Here 

the first term, the difference quotient 



I.e. 



A a a I stays finite in this limit. Further, 

dimension and indices are reduced. The remainder R4jj...k does no contain an inverse ()4. In a second 
step, in subsect. 5.2 we will eliminate the inverse ()4 in the difference quotient. 



5.1 Difference quotients with l/()4 

In this subsect. we derive optimized, compact expressions, where the appearance of possible singular 
1 / ()4-terms is reduced as much as possible. We will treat the singular behaviour using the fact that the 
integrals are exactly known in the hmit ()^ — 0. In fact, if ()^ = 0, due to ( |2.31| ) the n-point integrals 
degenerate to integrals with scratched propagators: 

n ft] 

Accordingly we define objects which converge in the limit ()4 — )• to the corresponding tensor 
coefficients, taken in that limit: 



5 (t^) 
f=l (0D5 



4'^" = ItSt^'^"' (5.3) 



- U lQs\ (Os\ M ' y-^-^) 

V 7W,-y _ V lo^^A/ Jd),st {oshjd),s {Qj_Jd),s .^rs 
f=Lf/;j- 1,0^5 Ws \0s)5 



5 (ts) 
\Os)t 

t=l,tj^i,j,k (0J5 
I (0^)5 ,, T{d),s,iLh rid),s iosh r{d),s 

+ 7o^^ijU,ij +7o^^ikh,ik +^Q^^jkh,jk- (5-6) 



V0J5 \Os)t 
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Eq. (|A.7|) reads in this notation 



j{d+2),s 
^4 



Os 
Os 

Os 
Os 



0. 

1 

(3; 



5 

E 



ios)5_r{d),st 



t=i ioDs 



d-3 



Jd),s y{d),s 



J d-3' 



such that indeed for — )■ the 74 remains finite. We need this relation with d 
d = [d+] = 6 — 2£ for tensors of rank i? = 2, 3 and rank = 3, 4, 5, respectively. 



(5.7) 



4 - 2£ and 



The next integral is The recursion of integrals with one index, 1)"/'% is (|A.6|). To rewrite <\AM 
in a similar manner as (5.7 1, we evaluate the right hand side of (A. 6), replacing I4 by Z4 : 



{d),s 



\isJ5 ' 







sJ5 



5 

r{d),s _ \isJ5 j{d),st 

'4 2^ (s\ ^3 

t=i U5 



1 



5 

f=i CDs 



1 



(0D5 ^= 



5 

0^/ 



(oJs 

lOiv' 5 



Jd),st 



Ad), St 



(5.8) 



where the latter eqn. is due to 

Os 
is 



Ost 
Osi 



(5.9) 



The factor 1 / has cancelled and we obtain the analogue to ( |5.7[ ): 



Jd+2),s 



©5 



Jd),s „{d),s 



+ ■ 



05? 
05/ 



{d-3)I, 



{d+2). 



+ E 

t=]. 



Ost 
Osi 



rid), St 
'3 



Ad), St 



(5.10) 



In the first 



ine of ( |5.10[ ) we have introduced a difference quotient which can in the next step be replaced, 
due to (15 .71), by 74^^^^ '\ i.e. an integral of the same dimension as the original integral on the left hand 



side. In fact the second line of (5.10) is already our final result for this type of tensor coefficient. We 
need this result for the 5-point tensors of rank R = 2 with d = 4 — 2£ (generic dimension), for 7? = 3, 4 
with d = [d+] , and for R = 5 with d = [d+] 
above. 



Eqn. ( |5.10[ ) demonstrates our principle as described 

In the tensor integrals of higher rank more complicated difference quotients appear and will be 
dealt with in the next subsect. |5.2[ The procedure of calculation is the same as before. In order to 



obtain, e.g., VijI^^^ '\ we calculate the right hand side of ( |A.5| ) (for / = 2), replacing I^- 
Using again relation (|5.9[) we find 



+lbyZl'^ 



4,/ 



i+],s 



'4,i 



^[d+],s 



+ 



V0J5 



Osi 
Osj 



j[d+],s 

I4 -r 



E 

t=l,t^i 



Ost 
Osj 



Ad+],st 



(5.11) 
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Next we obtain from (|A.4|) 



^ij^ijkU,iik 



\ks)5 

Os 

1 

\0s) 5 



"'J 



r[d+] 



-z 



(5.12) 



and with(A.3) 



(::)5 

0^/ 



yld+]\s' 



+ 



\OsJ t 



0^/ 



^ fOst 

t=l,t^i,j,k N / o 



^ij^ijkhjjk 



(5.13) 



We now have collected all contributions to higher-dimensional integrals with an (|^) ^ in the denominator 
in such a way that also the numerator vanishes for (^,) g = 0, see ( 5.2 1. These results are only a rewriting 



of the recurrence relations, but they make the finiteness of the integrals at (^) g = manifest. They will 
be a starting point to find a final representation, which is truly optimal for Idnematical points around 
()4 = 0. In the second line of (5.10) we observe that there are no explicit inverse Gram determinants 



anymore. In the following we will show that this also holds for integrals with any number of indices. 



5.2 Reduction of the difference quotients 



In app. |A] we reproduce a list of the recurrence relations needed for the evaluation of the 5-point 
functions. In fact, since all tensor coefficients of the 5-point functions have been reduced to higher- 
dimensional 4-point functions, we need only the recursions for the latter. When applying these formulae 



to 5-point functions, we have to identify ()^ 



Ogand/^ 



d+] 



, etc. In the present sect, we will 



drop the index s in the Gram determinant and in the upper indices of the integrals. 

We now discuss the higher-dimensional 4-point fu nction s needed for the different tensor r anks 



0). 



Ti 



■\d+] 



of the 5-point functions. For the tensor of rank R = 2 (4.21 ) we need ^ and ■ ^ given in (5.7 ) 
and (|57 

R = 3, 
and (|5 



n the spirit of our approach they are already in the final form. For the te nsor of rank 



(|4.36) and (|4.37 
In fact 



we further need /i 



I^f:^^^ and V;//r.T^ . These are given in (15. 7|), (15.1 01) 



4 ' 4,i " 'J 4Jj . . , 

the first two are already in the final form, while in the last one a new difference 
quotient appears. 

Our general approach to cancel Gram determinants is, first of all, to use the recurrence relations 
"backward", i.e. to express a 4-point function of dimension d by one of dimension d + 2, multiplied 
by a Gram determinant, plus a sum over 3 -point functions. The factorized Gram determinant can be 
cancelled and for the collected sum over 3-point functions the algebra of Cayley determinants allows 
to combine them such that again the Gram determinant factorizes and can be cancelled. 



With the notation () = ()4 we obtain 
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(5.14) 
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and from ( |5.11| ) 
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(5.15) 



This is the form we wanted to obtain. The higher-dimensional 3-point functions can be calculated 
by means of the recurrence relations given in app. |A.2[ reducing them to scalar functions in generic 
dimension. If () is not small, the same applies for the higher-dimensional 4-point functions, in 



Jd),t 



in (5.15 1 in the same 



particular ( |A.7[ ); otherwise we will use ( |5.3| ) by setting up an expansion in the small Gram determinant. 
This will be done in subsect. 15.31 

It is worth mentioning that we can deal with 3-point functions like I\" 
manner as we dealt with the 4-point functions: 
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(5.16) 



to be compared with ( |5.10[ ). This allows to handle the 3-point functions in case ([) = 0, for which 
case (A. 1 1 ) does not work - expanding in small (J) if needed by the use of (5.2 1. The (q^) , however, 
vanishes for an infrared 3-point function - thus we have to assume here that (^) and (q^) don't vanish 
simultaneously in order to be able to apply at least one of them; this is also implicitly assumed for the 
case of the 4-point functions. 

Exploiting this approach systematically, it can be achieved in general that the indices are carried, 
like in ( |5.15[ ) and ( |5.16[ ), only by the Cayley determinants, multiplied by scalar integrals in higher 
dimension. This property might become useful for further analytical evaluation of the original Feynman 
diagrams, in performing partial sums over indices explicitly where needed. We point out that due to 
the powers of d in front of higher-dimensional integrals we have to take into account finite rational 
contributions arising from the diverg encies o f the integ rals; s ee app. |B]for a list of examples. 

For the tensors of rank R = 4(see (4.45 1, (4.48 1 and (4.53 1) we further need 4"^^^ and Vij^ijkl^^^^^- ■ 
The tensor with two indices was treated m (5.15 1 and we have only to shift the dimension: d ^ d + 2. 

Much more involved is now the calculation of if^l . The crucial point of our approach is to obtain 
here an expression for the following difference quotient with the envisaged properties: 
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(5.17) 



By construction, the 4-point functions have no explic it inverse Gram determinant anymore. It is 
interesting that the formerly calculated if^} (see (5.15 1) enters here as a whole. The remaining task is 
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now to show that in the sum of 3-point functions in (5.17) a Gram determinant () factorizes and thus 
cancels its overall factor 1/ (). Indeed this is so. After a tremendeous amount of cancellations one gets 
the result: 
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(5.18) 



The vul^^ 

as the final result. Adding all contributions, using (|5.10[), we finally have 



from (5. 15 ) (J — 7- J + 2) has now been explicitly inserted since it has the same structure 
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(5.19) 



It is obvious where the various contributions come from: Those being proportional to — / (q) 
come from ( 5.18| ), while all the others come from the second part of (5.12). We could indeed list 
these contributions separately without inserting the second part explicitly. However, in this way the 

coefficient of I^^^ gets contributions from both terms which combine to make the resulting coefficient 
explicitly symmetric in all indices. The symmetry of the 3-point functions is not so easily seen. 
Neverheless, numerically, it might be faster not to combine these terms but just to retain what had been 
derived. 

Also here, as discussed in ( 5.16[ ), we can replace the tensor- 3-point functions: 
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(5.20) 



At this point, we observe a simple rule of how to obtain ( 5.19[ ): replace in ( 5.15| ) d ^ d + 2 and 



multiply with — (J — 1) where (d — l) is to be chosen such that all factors {d + i) increase by steps 

Ko) 



of 1 (see also ( 5.15^ ). In this manner we increase simulaneously the dimension and the number of 



indices. Then one has to add the second part of (5.12 ). 

r[d+]' 



Finally, for the tensor of rank i? = 5 we need 



4,i jkl ■ 



Due to the above rule we need not again 



perform a complicated calculation, rather we apply the rule as proceeding from (5.15) to (5.19), 



increasing simultaneously dimension and number of indices: we have to shift in(5.19)J— )'<i-f2. 
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(0) 

muliply with — (iv^Y ^^'^ ^^^^ ^he second part of (|5.13[). We obtain: 



0\ (^\ /'On /On 



(;•) ffl (?) + o (°) ffl + o © + (Si) (°) ffl + (S o 



d{d+\)r' 



(Oi) (Oj\ , (Oj)(0') I fO^UOi^ 
lO//* lO/t^ ^ lO/^ lOyt^ ^ lO//' lOj7 JrfH 



o^2 



(2) 



(°) ffl (?) rf(rf+ l)(rf + 2) ^ /'OA , (°) (?) + 1) f 

(1!)(1!)Q (S) kMO' ^00 h 



(S) a 
(2) (S) ^ 



0/ 
Ok 

Oj 
01 

Oj 
01 



Ot 
Oj 

Ot 
Oi 

Ot 
Ok 



+ 


k 



Oj 
Ok 

+ 



Ot 
Oj 







+ 



Ok 
01 



jm\t 



+ 



OA fOt 



01 



Ok 
01 



Ot 
Oj 



'3,i 



Ot 
Ok 



+ 



(0\ 

t=l,t^i,i 



j[d+]\t 



(5.21) 



Again we, first of all, mention that the 3-point function I^^l appearing here can as well be calculated 
like (|5.16[) and (|5.20[), essentially by taking over the 4-point result, shifting J — )■ J + 1 : 
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(5.22) 



In fact, having a closer look at our 4-point tensor coefficients, we observe that apart from the 
higher dimensional 4-point functions, all other terms are 3-point tensor coefficients, occasionally of 
higher tensor rank. Eq. 5.21| has been obtained by an educated guess. Indeed, a step by step derivation 
would have been extremely tedious if one would have had the courage at all to try the calculation. 
Of course one needs a verification by numerical checks: among others we found for non-exceptional 
Gram determinants an agreement with LoopTools of typically more than ten decimals; see sect. |C] 
for some details. Concerning these non-exceptional Gram determinants, we just mention that we 
evaluate ( |5.15p , ( |5.19p and ( |5.2ip by means of the recurrence relations of app. |Aj For further details 
see Sec.O 
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5.3 Expansion of I\ ^ for small Gram determinants 



The tensor coefficients in (2.171 to (2.20), in particular ( |5.15| ), ( 5.19| ) and (5.21 1, have been expressed 



in terms of 4-point functions in higher dimensions: ,1^ 



4 . In our approach 

only these integrals can cause problems for small Gram determinants and therefore, finding a special 
approach for their calculation, will finalize the problem of calculating the 4-poin tensor coefficients. 
We start from the fact that for exactly vanishing Gram determinants ( |5.3[ ) yields a finite value for 
4-point functions of any dimension and, taking into account higher orders, we set up an infinite series 
in terms of powers of the Gram determinant]^ 
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l,---,4, 



(5.23) 



where the coefficients /f ■ have to be determined and 



(5.24) 



We use the observation that the recurrence relation with shift of dimension ( |2.31[ ) (see also |A.7| ) for the 
scalar higher-dimensional 4-point functions can be written as follows: 



+ « [(2/+l)-2e]/ri""', 
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1 = 1. 



where we re- wrote ( |5.3[ ) as follows: 
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(5.26) 



From now on we assume Z > 0, so that the integrals are infrared finite and the leading singularity in e 
is at most of the order l/e. Integrals in generic dimension (/ = 0) will be discussed at the end of the 
section. We treat the finite and divergent parts separately: 
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The first few iterations of ( |5.25| ) are 



Ad^ 



+ 



I 4d 



= z 



= z 



'd+\ 



'd+\ 



+ " r 7.' 



(5.27) 
(5.28) 

(5.29) 



y[d+r' , ^ J, 
4 ~'" ToV '+2-'4 



0^ 



.\d+ 



1+2 







.ON -(+1-4 + —2C1+1C1+2 < Z4 + 7o\"^'+3^4 



= z 



- (2)' 



,7=1 



Ad- 



(5.30) 



■'The series is not a Taylor series because the expansion coefficients are not the derivatives of ^ at () — 0. 
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with 



ci+j = 2(Z + j)-l-2e. 



(5.31) 



The 4-point functions are expressed in terms of an infinite power series in () / (q) with higher- 
dimensional 3-point functions in the expansion coeffients. 



For the finite and divergent part of (5.27 ) we get to lowest order 
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(5.33) 



Now the second part in (5.32) is proportional to () and can be considered as a correction term for 

[rf I iC+i) 

small Gram determinants, where also a proper approximation for ^ has to be chosen. There 
are two simple choices. One has just to set in ( |5.32[ ) = Z^p and neglect the unknown higher 



We 



order terms, or one selects a close kinematical point with () = and uses = | q^q. 

come back to these two alternatives later. 

To calculate higher order corrections, we perform now iterations. Defining the correction term as 



(2/ + 1)2^+1) -24^+1^'"^ 
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(5.34) 



the iterative scheme then reads: 



(5.35) 



The index / counts the highest power of ()4 and the series ■ is supposed to converge for growing / 



towards F^ . As a condition of applicability of the iteration we can obviously use 
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where the scale has dimension of a squared mass. It is worth to perform the first few steps in the 
iteration explicitly: 
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Performing / steps in the iteration, we have 
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where r is given in (|5.24[) and 



r(L+i) 



(5.39) 



The T{z) is the Euler Gamma function. In (5.38) we have to define yet Z^^^'^. Strictly speaking for 



4i;^'^=Fr'"''wehaveZ« 



'4F,i 



(i= 1, 2 ■ ■ ■ ) but in order to evaluate (5.38 1 we have to choose 



an appropriate approximation. 

Taking / — )• oo, however, and assuming convergence of the series, we have 



(5.40) 



i.e. in this limit the term with Z^q^ drops out. The choice of an approximation for Z^^q'' thus can 
influence only the first few partial sums of ( |5.40[ ). In fact, as will be seen in Sec.|C]in an example, the 
convergence is quite good for moderate r and after a few steps the result is not very much dependent 
on the approximant of Z^q^ . 

It remains to deal with the last term in (5.38 1, i.e. £>4 .As in (5.38 1, we define a partial sum 
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and get corresponding to (5.40) 
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y I i(L+;+l) 

In order to avoid in ( 5.38 ) higher order terms coming from D)^ ^ , i.e. higher than those contained 
in the finite part, we re-write 



r ^A 



(5.43) 



E.g. for j = i — I on the right hand side of (5.43 ) there contributes a term r^Z^^^, where / is the 



highest power of r which occurs also in the finite part of ( |5.38| ). Some algebra yields 
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with coefficients 
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such that as final expression in terms of an infinite sum, we can write the solution of ( |5.23[ ): 



I 

7=0 



L = 0,--4. 



(5.46) 



The V/'(z) is the logarithmic derivative of the Gamma function (digamma function). 

It is interesting to note that ( |5.46[ ) is also valid for L = 0. P ossibl e add itional infrared divergent 
terms are then contained for j = in '^^\ Comparing (5.46) with (5.23), we see that our goal is 
achieved. 

Having now obtained the compact expression (5.46), it remains to discuss how to calculate the Z\. 



Quite naturally one first calculates the needed z]'^^^ by means of recurrence (|A.11|) and sums over 
t according to (5.26). This is possible since in general for ()4 = 0, {^^ ^ 



^ for f 

' ■'2 



0, 



,4. Thus in 



order to evaluate ( |5.46[ ), what is needed at the end, are the , ? = 1 , ■ • • , 4 , w = 1 , ■ ■ ■ , 4, ? 7^ m 
and 4 1 -p oint fu ncions for the kinematical point under consideration, i.e. we need 14 master integrals. 

Applying (A. 11 ) in order to get the finite parts of /g"^^' '\ one also has to calculate the divergent parts of 

These have been discussed in 
^ 10 and larger), and for larger 



'"and/f 



the higher-dimensional 2- and 3-point functions I2 
app. |B] but only for low values of Z. They are needed for quite large Z (Z 
Z the analytic expressions blow up considerably. Further, for Z > 6 the analytic cancellation of the 
occurring Gram determinants is hard to perform. So, we preferred to work numerically. Amazingly, the 
situation is different for the 2-point functi ons. W ithout any problem we can produce with Mathematica 
any higher divergences with recurrence (A. 12), cancelling thereby the (J^) Gram determinants. For 



details see the discussion in sect. A. 3 To remain numerically as accurate as possible we use these 
analytic expressions when calcula ting the divergences of the 3-point functions numerically by recursion. 



starting withD3^^'(?) = — 5, see (B.4). As mentioned above, this has to be done for Z = 1, 



As a result, in a practical calculation one has the objects 
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which are a sequence of approximations for the finite part of the integrals 
we give an example for their numerical evaluation. 



, I.e. 



At this point we stress that formulae (5.15), (5.19) and (5.21 ) are free of indexed inegrals I\ 

' * ' ' ' ' '4,/'' 



In app. 



c 



and thus enable a new access to the calculation of tensor 4-point functions. Relation ( |5.46[ ) was 
already obained in Q, see (36) there: for n-point functions {n = 1,2,3,4) of arbitrary dimension with 
generic indices V, = 1 this series was derived in similar manner as above and a general scheme was 
developed of how to find analyic continuations to kinematical domains, where this series does not 
converge. E.g., if (q)^ is small, r in ( |5.24[ ) is large. For this case [7] contains the description of how to 
modify the procedure of solving the recursion such that the expansion parameter is small, i.e. how to 
obtain from the recursion relation a series in ^, see eqn. (23) ibid. Beyond that, from this series the 
J-dimensional n-point functions are obtained iteratively in terms of multiple hypergeometric series 
with ratios of different signed minors as arguments. For the 4-point function, e.g., the generalized 
hypergeometric functions 2^1, Appell function Fi and the Lauricella-Saran funcion Fs appear, see 
e.g. (98) of Q. Transformation formulae of the generalized hypergeometric functions allow to extend 
their applicability to different domains of the phase space. In particular in our situation we have to deal 
with integrals of dimension d>6 — 2£. The hypergeometric functions in [7] are also expressed in terms 
of 1 -dimensional integrals and inspection shows that for the large dimensions these are particularly 
well suited for numerical evaluation (see eqns. (78) and (96)). We leave this for further study. 
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Another attempt to perform the described series of approximations was undertaken in [l32ll : see 
Eq. (5) there. A specific example was studied, namely forward light-by-light scattering through a 
massless fermion loop. The approach was then not further followed. 



6 Symmetrized recurrence relations 



So far in the former sect, we were concerned with the evaluation of the 4-point tensor coefficients 
for small Gram determinants. If, however, the Gram determinant is not small there are other ways of 
doing the reduction. In fact the "standard" Passarino-Veltman [|3J reduction is one possibility. This is, 
however, not a unique procedure. 



While in [21J a systematic application of recursion relations of type (2.5 1 was performed for all 
tensor n-point functions, here we take a different point of view, namely to arrange the tensor coefficients 



in (2.17 ) to (2.20) for the 4-point functions in such a way that a possible analytic simplification of the 
tensor as a whole is achieved. ^ 

We begin with ^ defined in ( |2.18[ ), containing as most complicated object if^^ . This integral, 
represented also by ( |5.15[ ), may be reduced by ( |2.30| ), 
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The first observation of interest here is the symmetry in the indices z, j. Only the last term is not 
obviously symmetric. As was mentioned earlier, the symmetry is in general seen only after summation 
over s,t, which we can exemplify here. We use the relation 
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(6.3) 



Inserting the left hand side into (6.2 1, the first term on the right hand side has just exchanged indices 
ij. In the second contribution cancels due to which the sum over s,t vanishes since i^J^ is 
antisy mme tric 'mt,u. 



In (16. 11) there remains as higher-dimensional integral lt^\ which is evaluated according to (|A.7|). 



Often it is as well used as "master integral" since it is UV and IR finite. Having a look at (2.18) we see 



that the second amphtude of the rank R = 2 tensor is alsqjust if^K This allows the following way of 



writing for this tensor. Similarly as in [|2T]| we introduce [j 
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''C^ here differs from the definition (24) in 121] by a factor of 2. 
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with 



(6.5) 



and = \{)- This allows to drop if^^ in (6.1 1 and to replace in (2.18 1 by G'^^: 
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This representation may become advantageous in the analytic evaluation of diagrams since G^^, 
contracted with a (proper difference of) chord(s), vanishes. Remember that any external momentum 
may be written as a sum of chords. For the tensor of rank R = 2 this corresponds to (20) of [1211 . 
We will derive here t he cor responding relations for the higher tensors as well. 

Proceeding to 7^''^, (2.19), we , first of all, need in the second term on the right hand side 1^^^'^ , 



((A.6) with I = 2): 



where again I^f^^ appears and /!^^^'' is given in (A.l 1 ). From the recursions of subsect. 
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we obtain 
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where we have introduced an abbreviation for the remaining 3- and 2-point functions: 
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In (6.8 1, the 







of (6.1 1 together with part of the last contribution in (A.9 1 has now been absorbed 
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: in (A.4) only two terms of this type appear explicitely, i.e. this form has no obvious 



symmetry. Further, in ( 6.8[ ) and ( 6.9^ only the 4- and 3- point functions are explicitly symmetric in 
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the indices ij^k. In order to demonstrate symmetry also for the 2-point functions one would have to 



reduce also /. 
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which is given in app. 



We may now combine the results and simplify I^^^ correspondingly: The /j"^^' in (6.8 1 can be 
combined with the second part of (2.19), i.e. it can be dropped in (6.8), and in (2.19) the must 
then be replac ed by G^^. (6.8 ) has also b een identified with a complete reduction of (5.19). 
For I^"^'' \2.20\, we start from (|aj|): 
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As can be seen, the ^ij^il^^^li can be expressed in a form which mainly contains terms which also 
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see (A. 5 ). The niji^I^^j. can be written similarly like (6.8 ): 
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where ^2^;^} collects the remaining 2- and 1-point functions and is obtained from (6.9 ) as follows: all 
3-point functions are replaced by 2-point functions (3 2) and the 2-point functions are replaced by 
1-point functions (2 — )■ 1). All summation indices u must be replaced by v and summation indices t 
must be replaced by u. Finally in all determinants and integrals columns, lines and propagators t must 
be scratched - like in (16. 121), i.e. 
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u,v=l \tJ \tuJ 



(6.13) 



and 



R 



[d+f,t 

2,ij 



(Ot) 4 fut\ 4 (ut\ 

(t\ 2^ (t\ n + 2^ (t\ h,i 



C) 



(6.14) 
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To indicate how the various tensor components may be combined to simplify the result, we only count 



here the number of contributions of a certain type. The second term in ( |2.20| ) contains, according 



to ( |2.23[ ), 6 terms of type ( |6.1 1[ ). Each type of the 3 terms in ( |6.1 1[ ) is also contained in the other tensor 



components in (2.20): 



■m 



-jyl^ ■ occurs 6 times in ( |6.11[ ), 



^4 



occurs 3 times in (16.1 11) and occurs 3 times in the 3'^^ term of (12.201) 



• Ef=i Tr4'^,-^' ' occurs 3 times in (|6.1 1]) and occurs 3 times in (|6.121). 



Thus, rewriting ( |6.4| ) as 



4 ('.) 







(6.15) 



one can convince oneself that g^^ cancels and it remains G'^^, which after contraction with a chord 
drops out. 

As we have seen, the above treatment of tensor coefficients requires at least one step of iteration of 
the recursion relation, because of which it is applicable only for non- vanishing Gram determinants. For 
these, however, it is very useful when we consider 4-point functions obtained by scratching one line of 
a 5-point function. In this case it is indeed possible to perform a cancellation of terms, which above 
still have the factor G^^. Also for this case, which we dealt with in sect. |3| the above presentation of 
the 4-point tensor coefficients can be applied - only scratching of one propagator has to be taken into 
account. 

Beyond that the present approach results in reductions which make the symmetry in the indices 
/, J, k, I more transparent and as a consequence yield certain blocks which can be calculated separately 
and combined to yield the complete tensor coefficients. 



7 Analytic simplifications for contractions of tensors with chords 

Before a numerical program for calculations is set up, it turns out to be advantageous to simplify 
Feynman diagrams analytically. A standard example is the following. If in the numerator of a 
Feynman integral a scalar product qt -kofa chord and an integration momentum occurs, this product is 
usually expressed in terms of the difference of two scalar propagators which can be cancelled against 
propagators in the denominator. Already in [21 J an alternative was indicated, making use of the fact 
that the contraction of a vector of the type ( |2.6| ) with a chord yields a simple expression: 



n-l (0) J 

?rGo=E^,-^;7r^ = -;^(i'/n-i'n.), /=!,.. .,n-l, (7.1) 

j=l \)n ^ 



and 

n-l 



('^■) 1 

qi ■ Gs = XI li^j^ = :^ (Sis - 5ns) : / = 1, . . . , n - 1, s=l,...n. (7.2) 

i=] y)n ^ 
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In ( |7.1[ ) and (7.2 1 ^„ = is assumed since only in this case 



(7.3) 



which is needed for their derivations. Thus, if the reduction relation (2.5 1 for /^i- is contracted 
with a qi^^, we can advantageously apply ( |7.1[ ) and (7.2 1 with n = 5. 

In case one considers a process with 5 external legs, one can choose from the very beginning <?5 = 
in the tensor integrals. If, however, the 5-point tensor is obtained by reducing a 6-point tensor, cases 
with <?5 7^ will occur. In order to be able to apply ( |7.1[ ) and (7.2), it is recommended to perform a 
shift of the integration momentum like k ^ k + q^, i.e. qi — > qi — q^. Such a shift is not a problem at 
aU, nevertheless it is interesting to see how this shift can be implemented in the formalism. The scalar 
integrals and the signed minors are invariant under the shift. We exemplify this for /j^^, writing 



'J=l 



with 



The sums 



+ 



+ ^^I^rv,,/rj^^-«I v,,75}', (7.4) 



!,!=1 



On 



(7.5) 



can be obtained by making use of 



j 

" r 



'j 



n 



ij 



ij 



" /n 



n 



!=1 



n,i ' 



0. 



(7.6) 
(7.7) 
(7.8) 



(7.9) 
(7.10) 



Thus one has 



',7=1 y=i 



n,ij 



i=l i=l 



(7.11) 



using the abbreviation q'^ = q^ — qH, and with 



n,i 

i=l 



1=1 



(7.12) 



38 



derived in the same manner, with ( |2.17p and ( |2.18| ) the standard result of shifting the integration 
momentum is obtained. The point is that the extra contributions obtained by the shift contain only 
integrals which were needed already in the unshifted integral so that the shift does not require the 
calculation of any new integral, see ( |7.5| ). 

Assume now again that we are dealing with the 5-point tensor and have q^^ = 0. The above trick 
to avoid an increase of the tensor rank can be applied to (14. 81) as well: contracting with ^ the first 
term yields the contribution qi ■ Qq given in ( |7.1| ), for the second term we have to find a formula for the 
scalar product qi ■ Q^. Indeed, 



4 





05 

(5,, - 55,) + 



{Yi5-Y55) 



(7.13) 



For the first term in (4.8 1 also another possibility exists, provided a contraction with a further vector is 
available. In such a case the first term on the right hand side of ( |4.1 1[ ), which shows up explicitly in the 
tensor components Eij,,,, yields a double-sum like 

E {qa ■ qi) {qb ■ qj) ^ = ' (o) 5 ^ ^ (^^5 - ^^55) (4. - ^5,) . (7.14) 



Further sums are obtained if the 4-point tensors are contracted. These are all of the type i-C- with 
line s scratched. We just list a few of them: 

4 /05^ 



£ qa ■ qj y . 



E iq^-qi)iqb-qj) 

and for 3-point functions 

4 ftS 



1 

~2 

2(5 

3 A 
2\s 



55 



5 

+ 2 



iYa5~Y55) 



1 



_ 1 
~ 2 



:qa ■ qb 



-Osi^ab^as + Sis 



(7.15) 
(7.16) 
(7.17) 
(7.18) 



[{l-das)5at-{l-d5s)55t] 



^ f ist\ 



2(5 
2(5 

St 
St 



St 
St 



l-4r)4.+ J (l-55,)55. 



5 

Os 
ts 



'l-55s)55t + 



(7.19) 
(7.20) 

(1-55^)55. [^(7.21) 
(7.22) 
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Even a quadrupel sum appears: 



jj*..)(...)C;)X;:)X;:)rK;).(:).(::). 

In fact, there are many more such sums. Our conclusion here is that in every scalar, which is obtained 
by contraction with chords the appearing sums can be evaluated analytically in order to yield compact 
expressions. This is due to the fact that the indices of the chords /, j,c ... are carried by signed minors 
while the integrals don't necessarily carry indices anymore. 



8 Conclusions 



We have developed a new approach to reduce tensorial one-loop n-point Feynman integrals based on an 
algebraic method elaborated in earlier papers. The approach is worked out up to 6-point tensors with 
rank R <6 and a rule is found how to extend the method to higher ranks. The first step was to reduce 
5 -point tensors up to rank 5 to 4-point tensor coefficients given in terms of higher-dimensional, indexed 

4-point functions. The latter are expressed in terms of higher-dimensional integrals I^^^'^ , L = 1 , ■ ■ ■ , 4, 



plus 3-point tensor coefficients in ( |5.10[ ), ( |5.15[ ), ( |5.19[ ) and ( |5.21| ). So far no Gram determinants ()4 



occur. Inverse powers of do occur if the integrals /j^^' are reduced to standard Aq,Bq,Co,Dq 



functions in generic dimension. For small this is avoided by using the expansion ( |5.46[ ) in positive 
powers of Application of Pade-approximants based on the e-algorithm to this expansion allows to 

calculate the I'f^^ in a simple manner to such a precision that the complete phase space is covered 
with high numerical precision. A numerical opensource code of the formulae derived in this article is 
under development. 



As a matter of fact, (5.46) is a special case of the general method developed in [VJ. Apart from this 



special series expansion, the integrals are expressed in in terms of multiple hypergeometric 
functions 2F1, Appell function Fi and Lauricella-Saran function Fs. In this context, it is a crucial 

property of eqns. ( |5.10[ ), ( |5.15| ), ( |5.19[ ) and ( |5.21[ ) to be free of integrals I^^} with indices larger than 
one. Using the notion of special functions allows a variety of options to adjust to various kinematical 
situations, and it might be interesting to explore their potential for a further improvement of numerical 
programs. 
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A Recursion relations for the reduction of higher-dimensional 4- 
point functions 



In this app. we quote explicitly the needed recurrence relations for the reduction of the 4-point functions. 
In [6J three recurrence relations are given: one reducing simultaneously index and dimension, one 
reducing only the dimension and another one reducing only the index. In the present work only the 
first two relations are used, i.e. ( |2.30[ ) and ( |2.31[ ). Special cases of ( |2.31[ ) are ( |A.7[ ), ( |A.11[ ) and ( |A.12[ ), 
i.e. these are the ones used to reduce the dimension of the scalar 4-, 3- and 2-point functions. All the 
other ones are special cases of ( |2.30[ ) and reduce the tensor indices. Reducing 4-point functions, 3- and 
2-point functions are generated, for which we also have to give the corresponding recurrence relations. 



A.l Reduction of 4-point integrals 

To learn about practical applications of the recurrence relations it is useful to investigate the most 
complicated one under consideration in some detail. Recall the definition given after (2.27 1, 
[d+] ^ = A + 2l-2£ = d + 2l. With ( [2301 ), 



^ijklU,ijkl 



Aj[d+?.t (/)4.[rf+P , (/)4 Jrf+f , OAj[d+? IX 



The integral contributes to the ranki? = 5 tensor coefficient Eiji^i„^, see ( |5.13[ ). Strictly speaking (A.l I 
is in this form valid only if all indices /, j, k are different. In case that some or all indices are equal 
there is no repetition of the same terms on the right hand side and the question is how to take into 
account this property in a general manner. Let us recall that in fact in the original integral (|2.21 1 we 

have to deal with > where = ^ij^ijk^ijki- Thus it is recommended to multiply (A.l I with 



VijVijk- Let us further introduce for the sum of the last three terms in (A.l I the notation [ijk]^ (with 
repetition) and [O^lllrf (without repetition). Then we have the following useful relation: 

vv -liikW)^^ {'l)4r[d+?,^ ii)4r[d+?,^ (z)4.[d+l' 

yijVijkiljkl^^j -[Ijkl^' + djk-^l^jl^ +Oik-^l^i,^ +Oij-^l^-- 



0. 



-Vjk—U,jk +^ik—\ik +^'J—\ij ■ 



0. 







0. 



4,0' 



(A.2) 



Thus, with riiji^ = VijVijk ( A.l[ ) reads 

j[d ' 



''ijklU,ijkl 



0. 



^ijkh,iik 



+ 



t=\,t^i,j.k 

Correspondingly we have 



()4 



[d+]\t 
ijk^3,ijk 



nijkl 



+ 



0. 



[d^ 



+ 



()4 



yikL 



, (/) 

AAk + 



0. 



.(A.3) 



^ijkh,ijk 







()4 
(") 



+ 



4 

E 



0. 



[dM' 

O4 



4 

3,; 



+ 



Qj-Ad- 
().^4,; 



O4 



I'-l 



+ ■ 



^ 

/=!,/// vJ4 



+ 1'-'/ 



+ 



©4, 



[d- 







(A.4) 



(A.5) 
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and directly from ( |2.30| ) and ( |2.31 



I we obtain: 



■\d+\ 











0. 



4 

''4 



^4 



/-I 



f=l 
4 

E 



0. 



4r[rf+ 

^3 



/-I 



1 



d + 2l-5 



(A.6) 
(A.7) 



In (A.7 1 we can put J = 4 for / = 1 since is UV and IR finite. Concerning the factors riij^i... and 



Vijki--- we see that the recursions work as if there were no such factors at all: each recursion eliminates 
one of the factors Vij.... This property continues to be valid for 3- and 2-point functions. 



A.2 Reduction of 3-point integrals 



riiiirL 



[d+]\t 



[d+]',t 



r[d+]',t 



M+tt 



(fa)4,, T[d+?,t , 



4 fut\ 

E UJ4 ., r[d+]^,tu 



{kt)4 T[d+]-,t I (fa)4 y[^+l 



C) 



C') 



04 

/On 

/Of 



ii=l,u^t,i 



C), 



+1'-'./ 



u=\,u^t (f) 



(if) 4 rl^H 



/OA 

O4 



d+ 



1-1 



Of special interest is the case / = 1: the /f is IR finite but UV infinite. The ^ is in any case UV 
finite, however, if it is IR divergent then the coefficient (*)^ is zero. Thus for the 3-point function we 

1/e, results in 



1 



J + 2/-4' 



, (A.8) 
(A.9) 
(A. 10) 
(A.ll) 



can put d = A. On the other hand, the UV divergence of I^^^^'' , coming from 12'ijy 



j[d+],t 

'3,UV 



- 1 / (2£) . Thus we have to keep the factor 1/ (J — 2) in (A. 1 1 1 when multipying 



A.3 Reduction of 2-point integrals 

For the 2-point functions, surprisingly enough, a number of peculiarities occur. Let us begin with the 
scalar integrals in arbitrary dimension. The recurrence relation quite generally reads 



(Qtu\ 

\Otuhj[dA 
(tu\ ■'2 
\tu) 4 



tu 



{Otu) 4 



H (tu\ 
= 1,V7^?,M \tu)4 



\d+ 



tuv 



1 



d + 2l-3 



(A.12) 



The (j")^ = —2{qi — qiiY = —2q^, ij' ^t,u (i' 7^ /) is independent of masses and is in particular 
the argument of the 2-point function, i.e. I2 = /2(mi,m2,<?^). Quite often — 0, which has to be 
considered separately. This situation occurs, e.g., in our case of calculating the higher-dimensional 
4-point functions where the corresponding 2-point functions are generated by the application of the 
recurrence relations for the 4-point functions. A more physical case, e.g., occurs if we consider 
radiation of a photon from an internal massive line. For q^ 7^ 0, nevertheless, we can apply (A. 12 1. In 
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that case it is first of all worth to investigate the 1-point functions. These depend only on one mass m 
and can be expressed for arbitrary dimension as 



(2m 



2^/ 



d{d + 2)---{d + 2l-2)' 



(A.13) 



with 



r(i-f) 

(m2)l-i 



m 



(A. 14) 



The question now arises how to implement these expressions in recurrence (A. 12 1. Here the following 
observation is useful: in ( A.12[ ), performing the sum over v, there are only contributions if all indices 
t, u, V are different. In the /i -functions the propagators with indices t, m, v are scratched. Let us call the 
remaining propagator w. Since all indices are r unnin g from 1 to 4 and all are different, theire sum is 



t — u — v and the mass in /i, (A.13 1, can be expressed as m = ^J^io-f-M-vjo-/ 



10. Thus w= 10 

where the Yij is defined in ( |2.10[ ). This property has in particular been used to calculate the divergent 

parts of the 2-point functions I^^^ for t,u= 1 , ■ ■ ■ , 4 as it was needed in Sec. 
Making use of 



5.3 



{mj - mj)I^'^^ {nii, niii ^q^ = 0)= I'f'^ (im) - if^ {niii 



(A.15) 



we obtain the following relation of special interest: 



^1 -hVj 



4" 



jtu 



(0) 



d-l 



(A.16) 



where q^ = {qt — qj')^. The second term is obviously UV finite and vanishing for equal masses. Also for 
q^ = this expression makes sense, the difference quotient resulting in the D50-function in LoopTools 
notation [25]. In general, however, and in particular for high dimensions one should start from ( 5.2[ ) 

i.e. we work in this case with 



in case of q = 0. For convenience we may drop indices t, u in I2 
2-point functions only: 



[d+],tu 



/f^' {mi,mf,q^ = 0) 



{mi,mf,q = 0) 



f=i [0)2 



(A. 17) 



and with q^ — 0: 



O2 



(-1)' 



2 2 ' 
mf — m f 



nii ^ mil 



(A.18) 



such that 



iA {mi,mii,q = Oj 



2 

mf ■ 



(nii) + 



m 



2M 



(A.19) 
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For the case m, = m,/ and = the ratio (0)2/ (0)2 ~ 0/0 does not yield correct results |^ and we have 
to choose a different approach, based on relation (29) of |l6l (see also H, IH ): 



+ ,(d+2) 



j{d) 



With 



the eqnn. (A. 20) yields 



rid) 
'1,1 



_r{d-2) 

Jd) 

2m2 1 • 



The latter eqn. is obtained from ( |A.14[ ), and finally 
For the vector integral the recursion reads 



j[d+]',tu 



({)tu\ 

KitulA j[d+]' ,tu 



nu\ '-2 



4 (vtu\ 
I KituJA j[d+Y'\tuv 

+ 2^ /tu\ h ' 
i'=l,v/f,M \tu) 4 



and similarly to (A. 16): 



j[d+],tii 
^24 



11 q^ 



(A.20) 



(A.21) 



(A.22) 



(A.23) 



(A.24) 



(A.25) 



which seems particularly useful for numerical evaluations even if q^ = 0, yielding the function 
DBO{0, m\, of LoopTools. To clarify the notation we mention that the integral occurs with 

the indices j = i and j = i' {i ^ i'), corresponding to the chords ^, and qif and the masses m, and m,/ of 
propagators and c,v, respectively!^ 
For nii = Mil = m and = we have 



I^j^^\m,m,q^ = 0) = —-I^'{m,m,q'^ = 0) 



d-2 
Am?- 



(A.26) 



Further we need the tensor coefficients 

(Qtu 



j[d+]',tu 



\jtu) 4 [d+]'-\tu \jtu) 4 [d+]'-\tuv \jtu) 4 [d+]'-\tu 7 1--, 

\tu) A v=l,vj^tM,i \tu) 4 \tu) 4 



(A.27) 



^In fact the result depends on the order of taking the hmits m, — m,/ ^ or ^ 0. 

'^For the comparison wih LoopTools. e.g., we put q2 = 0, m,/ = m^^ and qi = —p,mi = mf^ . 
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For all the occurring integrals on the right-hand side of ( A.27 1 are given above. For =0 and 

nii 7^ m,/ we now have 



= i^v-y'Kf + ^Mf(i' = 0) + ^Mfii' = 0), « / u. 

lo/2 lo/2 lo/2 

(A.28) 

Note that for / ^ j the first term does not exist and for / = j we have V,, = 2. Thus we have for / ^ j: 



'^^(^ = 0) = W'tfi^ = 0) + = 0), 

[0)2 {0)2 



and for / = j we have 



'if (<f = 0) = ^I'lP (n' = 0) + ^ :^/f (m,-), f ^ i 



(j) 



2j\dM^(2 



OJ2 ^ jd, 







02^ 



(A.29) 



(A.30) 



where due to (A.20) 



'1,11 



12 1 



'1: 



(A.31) 



with if given in (A. 14). To evaluate (A.29) and (A.30) we need 



dl 



with 



(A.32) 



dq^ 



■ = 0) = ^ {^2 - 2 K + m?) + 1 (m2 - m^) ^ ^ I (^2 ^ 0) ^ , (A.33) 



i.e. the second derivative -^^iq^ = 0) is needed in addition to the function DBO used also in 
LoopTools. 

Last but not least we ha ve to deal with q^ = and m, = m,/ = m, where we avoid the appearance of 
a ratio ^ as follows: (A.20) reads in this case 



'2i 



d-2 ,d 



Arn^ 



(A.34) 



the latter eqn. being obtained from (A. 26). Assuming ^ = 0, all integrals - are equal and we 
have Q 



[rf+]2 ^ l J-2 

3 4m2 
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/f(m). 



(A.35) 



B Divergent parts of higher-dimensional integrals 



In our final results ( |5.15 ), ( 5.19[ ) and (5.21 1 we obtained contributions of certain higher-dimensional 
integrals muliplied with polynomials in J = 4 — 2£ such that the 1/e parts of the UV divergent 4- 
and 3-point integrals combine with the e-powers of the polynomials to yield finite contributions. In a 
numerical approach these contributions have to be explicitly calculated, and for that purpose we list 
the infinite parts of those integrals and of scalar 2-point functions appearing in the reductions. As 
described before, for the calculation of corrections for small Gram determinants, we need 3-point 
integrals of even higher dimensions, which are, however, too complicated to be listed - apart from 



being difficult to obtain. They have been calculated iteratively in the numerics, see sect. 5.3 for details. 
The higher-dimensional integrals are in general UV divergent and we write them in the form 



(B.l) 



The terms D^n^ are obtained from the recurrence relations (see app. A), starting from lower dimensions. 
We mention that, whenever in the following list two of the occurring indices z, j, k, t, m, v are equal, the 
corresponding D's vanish. 
For 1 -point functions, it is 



D 



1 



2-4---2(/ + l) 



10— /— M— V, 10— f— w— V' 



(B.2) 



For i,j^t,u we have for 2-point functions: 
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■,t,u 
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1680 

+ ^Y^j{Yii + Yjj)+2Y^] 
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D 



ld+]\t„ 



\35Y:1 + SY^Yj, + 3YiYfj + SY^^f^ + 35Yf^ 
5Yij (77,3 ^ ^ly.^ ^ ^YuYl + lYl) 
67-5 {5Yi -f 47„-7,, + SY^) + 207,3 (7„ -f- 7,-,) + %Y^^ 



+ 
+ 



D 



'2.i 



^[37,, + 27,, + 7,,]. 



(B.3) 
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For /, j,kj^t we have for 3- and 4-point functions: 

r.ld+],t _ _1 

^3 - 2' 

^r''' = - 4?^ {3 + Yij + J-*) + Yjj [Y^ + J^;7 + Yj,) + 7,, (y,, + 7,, + 7,,)] 

[J^V {Yjj + + iY.k + Y,k)+ Yk, [Yu + Y,j) ] ] 

^3 

I^^'= (B.4) 



and 
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~24' 



D^t^ = 0, 



JrfH 



'4,;- 



1 

6' 



£>r' = -Y^[i'll+n2 + i'l3+n4 + i'22 + i23+i'24 + i'33+l'34 + K,4], 

^4^^^' = ^ {3 [7l 1 (Fi 1 + 712 + Fl3 + 1^14) + 1^22 (^2 + Y22 + Y23 +Y24) + 
Y33{Yl3 + Y23 + Y33 + Y34) + Y44{Yu + Y24 + Y34 + Y44)] + 

2 [Yn{Yi2 + Yn + Fu) + 713(1^13 + Y14 + Y23)+ 1^14(^4 + Y24 + Y34) + 

Y23 {Y12 + Y23 + Y24) + Y24{Y\2 + i^24 + ^34) + ^34(^^13 + ^23 + ^^34)] + 
Yx 1 (722 + 723 + 1^24 + 733 + 1^34) + 722 (7l3 + 7i4 + 733 + ^^34 + 744) + 
i33(7i2 + 7i4 + 724 + i44) + ^44(711 + 7i2 + 7i3 + 723) + 

^^12^^34 + ^^13^^24 + ^^14^^23} , 

" 24' 
l^lf = 0, 



rf^r = 0. (B.5) 



C A numerical example 

In order to demonstrate the use of our small Gram determinant expansion, we reproduce the numerics 



for the topology shown in Fig. C.l (b), which arises from the on-shell six -point topology of Fig. C.l 
(a). The example is taken from [2J, and first results obtained with our approach were reported in [.33.1 . 
In LoopTools [|25l conventions, the tensor coefficients are defined as follows: 

3 3 

^MvA = E i^/^if;vis:aAj7 + £(^,ivifa+^vA^//x+^lM^'v)^00/- (C.l) 

;j-,Z=l i=\ 
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For our conventions see ( |2.19[ ). The external momenta are assumed to be incoming: pi = Pg+,P2 ~ 
P^-:P3 = Pvii + Pu-,PA = Pe- + Pd- The inverse propagators are here Cj = [{k — qj)^ ^'^^ 
LoopTools conventions Cj = [{k + Kj^if- —mj]. The Ki are the internal momenta, expressible by 
the pr. Ki = puK2 = Ki+ pi.K^ = K2 + P3,K4 = 0. Then, with pj = 5/, (p/ + pj)^ = Sij, we set 
si2 = teii:S23 = Sfivu:S3 = Svu:S4 = tgj. The Corresponding tensor integrals are, in LoopTools ft25il 
notation: 



0,M|,0,0) 

The Gram determinant is: 

()4 = Z1(3) = 



(C.2) 



Det {2KiKj) 



and it vanishes if 



^ed ^ ^erf.crit 



Vu 



In terms of a dimensionless scaling parameter x, 

hd = {^~\~^)hd,cnU 

the Gram determinant becomes 

()4 = ~2 X Sfxvutefii^txvu ~ ^vw'rte^)- 



(C.3) 



(C.4) 



(C.5) 



(C.6) 



Following [I2I, we have chosen 



2 X 10'*GeV^ 
1 X 10^GeV^ 
-4 X lO^GeV^ 



(C.7) 



and get ted,cnt = — 6 x lO'^GeV^. For x=l, the Gram determinant becomes ()4 = —4.8 x 10^^ GeV^. 
We also need the modified Cayley determinant: 



2M| 



Ml 



Ml 


MI — Sjj^vii 


Ml 







Ml 







—ted 




—ted 






^ jlVutejl'^'^ ^ztelA "^^Vuted ''r S^yii{Syii tgj tg^) 
+ Mz [slu + (ted - t-e^f - 2Syu{ted + tg^l)] ■ 



(C.8) 



From ( |5.34| ) we see that a small-Gram determinant expansion will be applicable when the following 
dimensionless parameter becomes small: 



R 



i)4 

0. 



X 5, 



(C.9) 
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(a) 



(b) 



Figure C.l: (a) A six-point topology; (b) a four-point topology derived from (a). 



where 5 is a typical scale of the process. We choose here s — Suvu- Obviously, due to (C.6), x must 



also be small. Indeed, for e.g. x = 0.01 we have R = — 0.064j^ 



In Fig. C.2 we show the tensor coefficient Diiii(x), and in Tabs. C.l -C.2 the tensor coefficients 
Dim and Dm are tabulated for the region of interest of x. Because we assume in our formulae qi = 0, 
and in LoopTools ^4 = is assumed, one has to care about specific correspondences, in particular: 



^1111 = "2222 X4 2222' (^.10) 

Dm = n222 x4^2^- (C-11) 

The integral ^42222 contributes also to the rank R = 5 tensor coefficients E2222m, see ( |4.62| ), and /|^222 
to the rank R = 4 tensor coefficients £222/5 see (|4.48[). 



The numerics have been performed with Mathematica v.7.0, using LoopTools v.2.5 for the scalar 
4-, 3- and 2-point functions in generic dimension. Of course the first step to be done is to ver- 
ify ( |5.21[ ), ( |5.19[ ), and ( |5.15[ ) by comparing their numerical evaluation with results from LoopTools 
v.2.5 This verification is naturally only possible for non-vanishing Gram determinants, i.e. in 
this case all (higher-dimensional) integrals have been calculated by means of the recurrence relations 
given in app. [a] For x E (0.1, 1.0) we find in general an agreement to more than 10 decimals: in 

Tabs. 

In 



these cases are marked by /4 2222 ^^'^ ^^Aiii ' respectively. In the next step we have to 

is the highest available 



evaluate ( |5.2ip , ( |5.19[ ), and ( |5.15[ ) for small Gram determinants. If I 
index for which the Z]^p are calculated, we see that the upper limit / 



imax in ( |5.38| ) can be at most 
4 as reference 



imax = Imax — L. Sincc the integral of highest dimension is , we take the value L 
value. 

Tabs. C.l - |C.2 have been produced with a varying number of correction terms, specified by 
imax = Imax — 4, dependent on the size of Q^, specified in terms of the parameter x (see ( |C.6[ )). It is 
clear that for very small Gram determinants a few terms only have to be taken into account. Higher 
terms don't change anything beyond the shown accuracy. 

Due to our discussion following ( |5.47[ ) we have a certain choice for the value of Z^^^ for i = imax. 



see ( |5.38[ ). The choices are ()4 = and ()4 7^ 0, taken at the kinematics under consideration. Due 



to (C.6) the two options are parametrized here by x, specified in (C.5 ). In the tables we mark the 



^The dependence of R{x) on x is shown in Fig. 1 of ll34l . 

'"in fact in some cases we found a slightly better agreement with LoopTools v.2.5 than with LoopTools v.2.4. 
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Figure C.2: The tensor coefficient Dim^x); the data points for a least square fit are shown. Input: 
grey ( red) dots originate from the expansion at small x, black dots from the values at larger x. 



approximations correspondingly by and x. We specify the expansions as [exp p,imax], where Z^p 



L+i. 



max ) 

,0 



in (5.38) is calculated for p = and p = x, respectively. 



We have, however, still another option. With (5.47) we provide a sequence of partial sums S' of a 



series expansion for F^^^^ . This is exactly the input needed for the calculation of a Fade approximation 

for fJ^'^^ . We apply the £-algorithm for sequence transformations [[35l[36ll . It is described in detail 
in (371. The £-algorithm allows an efficient calculation of elements of the so-called £-table. The first 
column is zero, and the second one consists of the sequence S' = Z^^^, the convergence of which shall 
be improved. From the first two columns, the others are determined iteratively: 

£^_i\ = 0, (C.12) 

£q ^ — Z^pi, i = 0,- ■ ■ ,lmax~ L, (C.13) 



-^k 



The £-table and the Fade table are related: 



e(l = [k + i/k], (C.15) 

where the symbol [k + ij k\ stands for the degrees ^ + z of numerator and k of denominator polynomials 
of the corresponding Fade approximant \k + ij k\ . We took the choice k = lynax — L and: 

£(? ^ [k/k] ,,,,,. (C.16) 

In the tables we present the Fade approximants together with the corresponding sums, for x = as 
well as for jc 7^ 0: they are denoted by [pade 0,Zma.,:/2] and [pade x,«„,ax/2], respectively. In general 
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the Pade approximants provide a remarkable improvement of precision compared to the sums such 
that we even close up to the values provided by the non-small Gram determinant representation. It 
is also remarkable that there is only very little difference between the values obtained by iterations 
starting from jc = or from x^O. In fact the Pade starting with jc = is only slighly better than the 

one calculated with x j^O. The difference arises because the integrals f}^^^^ change much slower 
with X than the approximants Z^p ■ so that for a start atx = the latter are already much closer to the 
final values than those starting atx ^0. Nevertheless, after a few iterations the difference has already 

almost disappeared. At this point we want to remind the reader that we approximate the integrals I^^^ 
and they must indeed be very precise since there are considerable cancellations with the remaining 
contributions to the tensor coefficients under consideration. 

To discuss a few results given in Tabs. C.l| and C.2 we mention in particular the results for x = 0.05. 
In Tab. C.l we can assume the values of [pade 0,10] to be accurate to 10 decimals. One reason is 
that even Pade approximants up to [pade 0,13] (based on 27 terms in ( |5.38| ), not shown here) remain 
stable in the lO''' and 11^'^ decimal. With LoopTools we have an agreement of 8 decimals - thus it 
seems that this point is the value where the representation in terms of non-small Gram determinants 
starts to loose precision. We also see that there is even agreement beween [pade 0,10] and [pade x,10] 
up to 9 decimals. Just for curiosity we have also calculated the values for x = 0.1. Here LoopTools 
and also our calculation with the larger Gram determinants work perfectly so that we can say to be in 
the domain of larger Gram deteteminants. Nevertheless [pade 0,13] yields a precision of 3 decimals 
- while [exp 0,26] is off by 2 orders of magnitude. It is clear that for lower x we obtain good results 
with less terms in the expansion. The results in Tab. C.2 are of similar quality - in fact they are slightly 

better since the highest dimensional 4-point function is only I^^^ . 

Finally we show in Fig. C.2 the smooth behaviour of Diiii{x). The smoothness is in striking 
contrast to the complexity of its precise calculation. Indeed, in the critical domain the function is 
almost constant and one has to spend quite some effort to calculate it with high enough precision. 
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X 




3m Dim 


0. [exp 0,0] 


2.05969289730 E-10 


1.55594910118 E-10 


10-** [expx,2] 
[exp 0,2] 


2.05969289342 E-10 
2.05969289349 E-10 


1.55594909187 E-10 
1.55594909187 E-10 


10-"^ [expx,5] 
[exp 0,5] 


2.05965609497 E-10 
2.05965609495 E-10 


1.55585605343 E-10 
1.55585605343 E-10 


u.uui [exp U,DJ 
[exp x,6] 

M,2222 

£>iiii 


2.05932484381 E-10 


1.55501912433 E-10 


2.02292295240 E-10 
2.01707671668 E-10 


1.54974785467 E-10 
1.62587142251 E-10 


0.005 [exp 0,6] 
[pade 0,3] 
[exp x,6] 
[pade x,3] 

M,2222 

£>iiii 


2.05786054801 E-10 
2.05785198947 E-10 


1.55131031024 E-10 
1.55131031003 E-10 


2.05786364440 E-10 
2.05785199805 E-10 


1.55131031024 E-10 
1.55131030706 E-10 


2.057788941 14 E-10 
2.05779811490 E-10 


1.55135794453 E-10 
1.55136343923 E-10 


0.01 [exp 0,6] 
[pade 0,3] 
[exp 0,10] 
[pade 0,5] 
[exp x,10] 
[pade x,5] 

M,2222 

£>iiii 


2.05703298143 E-10 
2.05600940065 E-10 


1.54669910676 E-10 
1.54669907784 E-10 


2.05600964693 E-10 
2.05600955381 E-10 


1.54669910676 E-10 
1. 546699 10676E- 10 


2.05600963675 E-10 
2.05600955381 E-10 


1.54669910676 E-10 
1.54669910676 E-10 


2.05600013702 E-10 
2.05600239280 E-10 


1.54670651917 E-10 
1.54670771210 E-10 


0.05 [exp 0,6] 
[pade 0,3] 
[exp 0,20] 
[pade 0,10] 

1 o V n V /III 

[exp x,zuj 
[pade x,10] 

j[d+]' 
''4,2222 

£>iiii 


4.83822963052 E-09 
2.01518061131 E-10 


1.51077429118 E-10 
1.50591643209 E-10 


2.04218962072 E-10 
2.04122727654 E-10 


1.51077424143 E-10 
1.51077424149 E-10 


2.04122727971 E-10 


1 S 1 077zl9zl 1 /l^ P 10 

1.51077423985 E-10 


2.04122726387 E-10 
2.04122726601 E-10 


1.51077422901 E-10 
1.51077423320 E-10 


0.1 [exp 0,26] 
[pade 0,13] 
[exp x,26] 
[pade x,13] 

'4,2222 

Dnn 


2.20215264409 E-08 
2.01749674352 E-10 


1.46815247004 E-10 
1.46681287362 E-10 


2.08190721550 E-08 
2.03995221326 E-10 


1.46815247004 E-10 
1.46785977364 E-10 


2.02269485177 E-10 
2.02269485217 E-10 


1.46815247061 E-10 
1.46815247051 E-10 


^ ■ '4,2222 

£>iiii 


1.72115440143 E-10 
1.72115440148 E-10 


9.74550747662 E-11 
9.74550747662 E-11 



Table C.l: Numerical values for the tensor coefficient Dim. Values marked by Dim evaluated 



with LoopTools, the I\ 2222 corresponds to ( |5.21[ ). The labels [ exp 0,2n] and [ pade 0,n] denote iteration 
2n and Pade approximant [n, n] when the small Gram determinant expansion starts at x = 0, and [exp 
x,2n] and [pade x,n] are the corresponding numbers for an expansion starting at x. 
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X 




3m Dill 


[exp 0,0] 


-3.15407250453 E-10 


-3.31837792634 E-10 


10-** [expx,l] 
[exp 0,1] 


-3.15407250057 E-10 
-3.15407250057 E-10 


-3.31837790700 E-10 
-3.31837790700 E-10 


10-"^ [expx,4] 
[exp 0,4] 


-3.15403282194 E-10 
-3.15403282194 E-10 


-3.31818461838 E-10 
-3.31818461838 E-10 


n nni \f^^r\ y fA 
u.uui L*^^P 

[exp 0,6] 

j\d+? 
M.222 

£>iii 


'K 1 S^fi7SzlSzl9Q V in 
—J. 1 J^J^Z.y J_.-iW 

-3.15367545429 E-10 


— J.J LOH^Jo / i.J\J lL- LyJ 

-3.31644587150 E-10 


-3.15372092999 E-10 


-3.31645245644 E-10 


-3.15372823537 E-10 


-3.31635736868 E-10 


0.005 [exp x,6] 
[pade x,3] 
[exp 0,6] 
[pade 0,3] 

M.222 

£>iii 


-3.15208222856 E-10 
-3.15208230282 E-10 


-3.30874035862 E-10 
-3.30874035931 E-10 


-3.15208224867 E-10 
-3.15208230411 E-10 


-3.30874035862 E-10 
-3.30874035867 E-10 


-3.15208269791 E-10 
-3.15208264077 E-10 


-3.308740061 10 E-10 
-3.30874002667 E-10 


0.01 [exp 0,6] 
[pade 0,3] 
[exp 0,10] 
[pade 0,5] 
[exp x,10] 
[pade x,5] 

M.222 

£>iii 


-3.15006665284 E-10 
-3.15007977830 E-10 


-3.299159261 10 E-10 
-3.29915888075 E-10 


-3.15007991203 E-10 
-3.15007991324 E-10 


-3.299159261 10 E-10 
-3.299159261 10 E-10 


-3.15007991217 E-10 
-3.15007991324 E-10 


-3.299159261 10 E-10 
-3.299159361 10 E-10 


-3.15008000292 E-10 
-3.15008000292 E-10 


-3.29915916848 E-10 
-3.29915915368 E-10 


0.05 [exp 0,6] 
[pade 0,3] 
[exp 0,20] 
[pade 0,10] 

1 V n v /I 1 1 
[cXp A,ZUJ 

[pade x,10] 

j[d+? 
4.222 

£>iii 


-1.34278470211 E-11 
-3.13432516570 E-10 


-3.22448580722 E-10 
-3.22580791799 E-10 


-3.13359445767 E-10 
-3.13365675001 E-10 


-3.22448581032 E-10 
-3.22448581024 E-10 


1 '^'^^T^n97izi p in 

— 0. V J\JZ.Z. LH- lL- LyJ 

-3.13365674956 E-10 


-3.22448581051 E-10 


-3.13365675084 E-10 
-3.13365675070 E-10 


-3.22448581110 E-10 
-3.22448581084 E-10 


0.1 [exp 0,26] 
[pade 0,13] 
[exp x,26] 
[pade x,13] 

j[d+? 
M.222 

£'111 


-2.49466252165 E-09 
-3.11144777695 E-10 


-3.13582331984 E-10 
-3.13599283949 E-10 


-2.34010823441 E-09 
-3.10806582023 E-10 


-3.135823319836 E-10 
-3.135870111996 E-10 


-3.11226750699 E-10 
-3.11226750695 E-10 


-3.13582331977 E-10 
-3.13582331978 E-10 


^- M.222 

£>iii 


-2.70193791372 E-10 
-2.70193791373 E-10 


-2.10251973821 E-10 
-2.10251973821 E-10 



Table C.2: Numerical values for the tensor coefficient Dm. Values marked by Dm are evaluated with 



LoopTools, the 222 ''^ defined in ( |5.19[ ). The labels [exp 0,2n] and [pade 0,n] denote iteration In and 
Pade approximant when the small Gram determinant expansion starts at x = 0, and [exp x,2n] 
and [pade x,n] are the corresponding numbers for an expansion starting at x. 
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D Notations and algebraic relations 

In the following we quote some relations for the signed minors, which are of particular relevance for 
the present work (see [|19.I and also IfTTll ). First of all we used the sums 

t(^) =0n (D.l) 



(=1 \'/ n 



and 



n 



i=l \'/ n 



0, (jVO). (D.2) 



Next a kind of standard relation, which allows to disentangle signed minors: 

and finally, what we call "master formula" in the present context ( (A. 13) of [[T9ll ) 

:)io:)„Ka,(::)/(:)„C:)-— ■ 

The following correspondences have to be taken into account when comparing our notations to that 

of e.g. m- 

()5 = -X^\ (D.5) 

= detX^^l (D.6) 



0/5 
0^ 



i/5 



4i\ (D-7) 



_ v(4) 

07 



X^j\ (D.9) 



- - ^^^^ mw^ 



^kjj^ (OO(^) Uk){Oi 

The notation in terms of D... -functions for 4-point tensors with even number of indices reads 



= l^kn^JJcj' =--- + g^'Doo. (D.ll) 
^n"^' = /£^^^^'^^'^^n^7' = ---+^"^V'''^oooo etc., (D.12) 



i.e. the functions Dqq... correspond in our notation to integrals in higher dimension /]^^^ (/ = 1, 2, ■ ■ ■ ) 
and our Z^f^^ -functions correspond to approximations of the Z)oo - -functions. We write 

Z^'+]' = zi^ + ^ + ff{£^) (D.13) 



and the series ( |5.46[ ) thus uses only approximations for the Dqo- -functions. 
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